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SECTION  1 


INTRODUCTION 


This  Research  and  Development  Technical  report  covers  the  work 
performed  from  1 June  to  31  December  1976  under  Contract  F33615-76-R-1216. 
This  report  is  submitted  in  fulfillment  of  CDRL  Sequence  3,  dated 
3 December  1975. 

In  accordance  with  the  Statement  of  Work,  the  overall  objective 
of  this  effort  is  to  develop  the  tools  and  techniques  required  to  assess 
the  technical  feasibility  of  using  a minimum  number  of  inertial  sensors 
to  provide  the  inertial-reference  information  consistent  with  the  re- 
quirements for  weapon  delivery,  flight  control,  navigation,  fire  con- 
trol, and  flight  safety  of  an  aircraft.  The  term  Multifunction  Inertial 
Reference  Assembly  (MIRA)  will  be  used  to  describe  the  generic  class  of 
systems  designed  to  meet  this  objective. 

This  MIRA  effort  is  primarily  oriented  toward  developing  and 
demonstrating  simulation  capabilities  to  evalute  MIRA-type  systems 
in  realistic  environments,  and  toward  consideration  of  the  life-cycle 
cost  implications  of  the  MIRA  approach.  The  attempt  is  not  to  propose 
or  evaluate  specific  systems  as  candidate  MIRA  approaches,  but  to 
develop  a versatile  evaluation  capability  which  will  have  long-term 
utility  in  the  MIRA  effort. 

The  effort  under  Contract  F33615-76-R-1216  is  divided  into  three 
tasks:  definition  of  program  specifics,  analysis  and  simulation  develop- 
ment, and  simulation  and  mechanization  studies.  The  results  of  the 
first  task  are  discussed  in  Section  2 of  this  report.  The  results  of 
the  second  task  are  discussed  in  Section  3 (for  the  aircraft  simulation), 
and  Section  5 (for  life-cycle  cost  models).  The  results  of  the  third 
task  are  discussed  in  Section  4 (for  the  MIRA  engineering  and  perform- 
ance studies) . Section  6 summarizes  the  results  of  the  effort,  and 
provides  transitional  information  related  to  utilization  in  the  develop- 
ment of  the  MIRA  concept. 


1 


SECTION  2 


DEFINITION  OF  AIRCRAFT  AND  FLIGHT-CONTROL  REQUIREMENTS 


2 . 1 Introduction 

At  the  beginning  of  the  MIRA  effort,  a survey  of  potential  candi- 
date vehicles  for  use  in  simulating  the  "high-performance"  application 
of  the  MIRA  system  was  conducted.  Although  it  was  recognized  at  the 
time  that  a model  of  a state-of-the-art  aircraft  such  as  the  F-15 
would  have  been  desirable,  it  was  felt  that  the  uncertainty  associated 
with  obtaining  sufficient  data  on  such  an  aircraft  was  large  enough  to 
seriously  jeopardize  this  relatively  short-term  effort. 

The  Charles  Stark  Draper  Laboratory,  Inc.  (CSDL) , with  the  coop- 
eration of  the  Air  Force,  obtained  both  rigid-body  and  bending-mode 
data  for  a typical  high-performance  vehicle.  The  data  were  sufficient 
to  allow  the  detailed  modeling  of  the  aircraft  dynamics  and  bending 
modes  (up  to  approximately  141  radians/second)  in  a real-time  all- 
digital  simulation.  The  availability  of  the  bending-mode  data  is 
particularly  important  since  it  is  anticipated  that  the  aircraft  system 
responses  to  the  bending  dynamics  will  be  a primary  contributor  to  the 
MIRA-system  location  sensitivity. 

The  overall  MIRA  program  considers  both  fighter-  and  transport- 
type  aircraft.  The  fighter  application  produces  a more  severe  environ- 
ment in  terms  of  accelerations,  angular  rates,  and  vibrational  excitation. 
Therefore,  the  limitation  of  this  MIRA  effort  to  a fighter-type  aircraft 
does  not  significantly  limit  the  utility  of  the  results  or  the  validity 
of  the  demonstration  of  the  simulation  capabilities. 

The  entire  simulation  is  structured  for  ease  of  modification  of 
the  aircraft  characteristics.  It  is  felt  that  this  capability  is  an 
important  feature  for  the  MIRA  program,  not  only  because  of  the  general 
nature  of  the  effort,  but  also  because  of  the  early  point  of  maturity 
of  the  program.  While  the  existing  simulation  utilizes  a significant 
portion  of  the  available  computational  resources,  it  is  felt  that  simu- 
lations appropriate  for  state-of-the-art  fighter  or  transport  aircraft 
could  be  developed  within  the  existing  simulation  structure. 
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2.2 


Aircraft  Definition 


The  data  selected  to  define  the  vehicle  were  the  most  comprehen- 
sive that  could  be  found  so  that  the  simulation  developed  would  have  a 
great  deal  of  flexibility  and  applicability.  For  example,  the  rigid- 
body  aerodynamic  data  defining  the  vehicle  were  developed  from  wind- 
tunnel  data  during  a program  investigating  high  angle-of-attack  flight. 
The  aerodynamics  are  available  as  nonlinear  polynomial  functions  valid 
beyond  the  stall  angle  of  attack  and  over  a wide  range  of  sideslip 
angles.  This  is  well  beyond  the  flight  envelope  needed  for  the  MIRA 
program.  The  data  are  valid  for  a vehicle  in  the  clean  configuration 
in  subsonic  flight  and  available  in  Reference  1. 

The  aeroelastic  data  defining  the  structural  modes  and  their  ef- 
fect on  the  vehicle  sensor  responses  were  obtained  from  References  2 
and  3.  Definition  of  the  vehicle  in  the  MIRA  simulation  includes  three 
longitudinal  and  three  lateral  directional  modes.  Each  mode  is  repre- 
sented by  a coupled  linear  second-order  differential  equation.  The 
longitudinal  modes  of  interest  are  the  stabilator-bending  mode,  the 
first  vertical-bending  mode,  and  the  stabi lator-rotation  mode.  These 
modes  were  selected  in  References  2 and  3 because  they  involve  strong 
f uselage-stabilator  coupling.  They  are  excited  by  the  stabilator- 
control-system  inputs  and  contribute  significantly  to  the  sensor  respon- 
ses. The  nominal  frequency  of  the  highest  mode  is  141  radians/second. 
The  lateral-directional  modes  included  in  the  MIRA  simulation  are  the 
fuselage  first  torsion  mode,  the  wing  first  asymetric-bending  mode, 
and  the  fuselage  first  lateral-bending  mode.  The  highest  mode  frequency 
is  80.8  radians/second.  Also  incorporated  into  the  MIRA  simulation  is 
the  effect  of  the  aeroelastic  modes  on  the  sensor  responses. 

2 . 3 Flight-Control  Definition 

The  flight-control  system  selected  for  primary  consideration  dur- 
ing the  development  phase  of  the  MIRA  simulation  was  designed  for  the 
Survivable  Flight-Control  System  (SFCS)  Program  and  is  presented  in 
Reference  2.  This  flight-control  system  was  developed  using  the 
structural-mode  data  programmed  in  the  simulation.  Together,  these  tv/o 
sets  of  data  provide  a consistent  package.  However,  the  flight-control 
system  developed  for  the  SFCS  Program  is  much  more  sophisticated  than 
needed  for  MIRA.  It  is  an  adaptive  system  valid  over  the  whole  flight 
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envelope  of  the  vehicle.  For  the  MIRA  effort,  the  adaptive  system 
was  reduced  to  a fixed-gain  system  suitable  for  operation  in  the 
vicinity  of  a nominal  flight  condition. 

The  flight-control  system  incorporates  conventional  feedback 
loops  to  improve  the  handling  qualities  of  the  bare  airframe.  For 
example,  the  longitudinal  system  incorporates  pitch-rate  and  normal- 
acceleration  feedback  to  the  stabilator.  The  lateral-directional  sys- 
tem incorporates  roll-rate  feedback  to  the  ailerons  and  yaw-rate  and 
lateral-acceleration  feedback  to  the  rudder.  A roll-to-yaw  crossfeed 
network  is  also  included  for  turn  coordination. 

During  the  course  of  the  simulation  dev'elopment , a second  longi- 
tudinal flight-control  system  was  programmed  into  it.  The  system  was 
obtained  from  Reference  4.  ilowever,  it  was  not  given  major  consid- 
eration during  this  phase  of  the  MIRA  effort  since  it  was  not  well 
defined  nor  consistent  with  the  aeroelastic  data  used,  as  in  the  case 
of  the  previously  discussed  system.  However,  this  system  is  avail- 
able in  the  non-real-time  simulation  should  it  be  desired  to  exercise 
an  alternate  configuration. 

2 . 4 Future  Aircraft  and  Flight-Control  Capabilities 

As  previously  mentioned,  the  basic  simulation  developed  for  the 
MIRA  program  is  capable  of  accommodating  a wide  range  of  aircraft  and 
flight-control  configurations.  Although  the  impact  of  the  MIRA  system 
on  the  flight-control  system  is  not  of  primary  concern,  and  a relatively 
straightforward  stability-augmentation  control  system  was  developed 
for  the  first  phase,  the  simulation  is  capable  of  supporting  consider- 
able sophistication  in  the  flight-control  area.  Among  potential  flight- 
control  developments  consistent  with  the  MIRA  approach  are  control  sys- 
tems for  Control  Configured  Vehicles  (CCV) , distributed  and  fault- 
tolerant  control  systems,  and  active  bending-mode  compensation.  The 
simulation  is  capable  of  supporting  these  efforts  so  long  as  the  control 
rates  are  maintained  within  a realistic  range  for  avionics  computers  in 
the  early  1980  time  period. 

In  addition  to  the  simulation  developed  under  this  effort, 

CSDL  has  a substantial  amount  of  data  on  the  KC-135  and  F-8  aircraft. 
Although  the  F-8  may  not  be  considered  an  appropriate  aircraft  for  the 
MIEIA  effort,  the  KC-135  could  serve  as  an  appropriate  model  for  the 
transport  application  of  MIRA,  thus  minimizing  the  program  dependence 
upon  external  data  sources.  The  reduced  dynamic  environment  (linear 
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and  angular  accelerations)  of  a transport  application  would  allow 
reduction  of  the  flight-control  and  environment- integration  rates,  and 
therefore  a reduction  in  the  overall  simulator  duty  cycle. 

Preliminary  investigations  were  made  during  the  course  of  this 
effort  into  obtaining  the  appropriate  F-15  aircraft  modeling  data 
for  upgrading  the  vehicle  model  to  a state-of-the-art  aircraft.  At 
the  current  time,  it  appears  that  sufficient  information  is  available 
from  McDonnell  Douglas  Aircraft  Corporation  (McDAC)  to  allow  modeling 
of  the  rigid-body  dynamics  of  this  aircraft,  but  not  the  flexible- 
body  dynamics.  Although  there  may  be  limited  benefits  in  doing  so, 
it  would  be  possible  to  convert  the  rigid-body  dynamics  to  an  F-15 
while  using  nominal  (or  possibly  new)  data  for  bending  dynamics.  This 
approach  might  allow  closer  correlation  of  the  CSDL  effort  to  the 
McDAC  MIRA  effort;  however,  it  represents  a relatively  inefficient 
utilization  of  the  simulation  since  the  bending  dynamics  represent  the 
highest  rate  computational  load,  and  can  have  a significant  impact  upon 
the  MIRA-system  performance.  Presently,  the  simulation  developed  under 
this  contract  can  run  in  real  time  without  the  bending  dynamics.  A 
real-time  simulation  with  bending  dynamics  will  be  available  early  in 
1977  in  the  second  phase  of  the  MIRA  effort. 

Although  the  choice  of  aircraft  model  involves  a wide  variety  of 
possible  options,  the  utility  of  the  existing  model  should  not  be  un- 
derestimated. The  dynamics  are  representative  of  a large  class  of 
fighter  aircraft  which  are  of  interest,  while  exact  aircraft  char- 
acteristics are  not  very  significant  this  early  in  the  MIRA  program. 

In  addition,  there  is  a large  body  of  information  available  concerning 
the  vehicle  which  may  be  useful  for  comparative  tradeoff  studies  of 
potential  MIRA  configurations. 

' 2.5  Inertial-Sensor  Requirements  for  MIRA 

2.5.1  Introduction 

Perhaps  the  most  disparate  requirements  for  inertial  sensors  in 
a tactical  aircraft  are  those  imposed  by  inertial  navigation  on  the  one 
hand,  and  stability  augmentation  on  the  other.  In  the  former  case, 
long-term  accuracy  of  the  inertial  components  per  se  is  the  prime 
requisite.  The  data  from  the  inertial  sensors  are  digital  and  repre- 
sent the  integrals  of  both  specific  force  and  angular  rate  (in  a strap- 
down  inertial  measuring  unit  (IMU)).  In  the  latter  case,  accuracy 
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requirements  on  the  inertial  sensors  are  orders  of  magnitude  less 
stringent.  The  data  from  the  sensors  are  analog  and  represent  the 
lateral  and  normal  components  of  the  specific  force  and  angular  rate 
about  the  roll,  pitch,  and  yaw  axes  of  the  vehicle.  In  the  absence 
of  a separate  inertial  navigation  system  (INS),  attitude  and  heading 
reference  data  are  supplied  by  a two-axis  free  vertical  gyro,  and  head- 
ing by  a single-axis  directional  gyro. 

For  further  tactical  applications,  the  prime  stability-augmentation 
sensors  (i.e.,  the  two  accelerometers  and  three  rate  gyros)  will  be  at 
least  triply  redundant.  Hence,  by  the  addition  of  a longitudinal 
accelerometer  to  each  of  the  three  sets  of  sensors,  one  has  all  the 
sensors  requisite  to  a triply  redundant  strapdown  IMU.  Of  course,  the 
stability  and  accuracy  of  the  individual  sensors  would  have  to  be  up- 
graded to  be  compatible  with  the  inertial-navigation  requirements. 

The  resulting  configuration  consisting  of  three-gyro-plus-three- 
accelerometer  prealigned  inertial-sensor  assemblies  could  be  distributed 
throughout  the  aircraft  in  positions  approximating  the  optimum  locations 
for  the  individual  stability-augmentation  sensors.  It  is  obvious  that 
this  configuration  contains  the  minimum  number  of  components  for  a given 
degree  of  redundancy,  and  hence  should  be  the  prime-candidate  configura- 
tion. 

2.5.2  Approach  for  Determining  Sensor  Characteristics 

Determination  of  the  characteristics  of  the  inertial  sensors 
which  simultaneously  satisfy  the  requirements  of  the  stability- 
augmentation  and  the  inertial-navigation  systems  necessitates  a multi- 
stage approach.  First,  the  stability-augmentation  or  flight-control- 
system  sensor  requirements  are  established  from  considerations  of  the 
vehicle  dynamics,  e.g.,  maximum  roll  rate  and  normal  acceleration,  the 
natural  frequency  and  damping  ratio  of  the  typical  (analog)  rate-gyro 
and  accelerometer  loops,  and  tlie  threshold  and  linearity  of  the  outputs. 
It  may  then  be  necessary  to  reinterpret  the  servo  characteristics  in 
the  time  domain  as  time  delays.  Next,  the  INS  requirements  imposed  by 
target  acquisition,  fire  control  and/or  weapon  delivery  are  considered, 
and  the  necessary  position,  velocity,  attitude  and  heading,  and  angular- 
rate  accuracies  are  determined.  Other  factors  include  reaction  time, 
radio-navigation  aids,  and  air-data  system.  The  range  of  INS  parameter 
accuracies  is  tabulated  and  typical  requirements  are  formulated.  The 
INS  requirements  are  interpreted  in  terms  of  inertial-component 
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performance,  stability-augmentation  requirements  are  superimposed,  and 
the  literature  is  searched  for  systems  and/or  components  which  satisfy 
these  composite  requirements.  Finally,  dynamic  and  performance  models 
of  the  inertial  components  which  satisfy  the  composite  requirements  are 
found  or  constructed  for  use  in  the  system  simulation.  It  should  be 
acknowledged  at  the  outset  that  it  is  impractical  to  consider  mechani- 
zing the  integration  of  the  exact  equations  of  motion  of  the  accelerom- 
eter pendulum,  which  has  a time  constant  ranging  from  a few  microseconds 
to  a few  hundred  microseconds,  while  the  natural  frequency  of  the 
stability-augmentation-system  outer  loops  is  on  the  order  of  1 hertz, 
and  the  natural  frequency  of  the  INS  is  0.0002  hertz.  Interest  in  the 
dynamic  sensor  performance  extends  only  to  frequencies  on  the  order  of 
20  hertz. 

2.5.3  Inertial-Sensor  Requirements  for  Stability  Augmentation 

The  stability-augmentation  sensors  are  the  normal  and  lateral 
accelerometers  and  the  roll,  pitch,  and  yaw  rate  gyros.  These  instru- 
ments provide  analog  voltage  outputs  proportional  to  specific  force 
and  angular  rate  respectively.  Restoring  force  (torque)  may  be  pro- 
vided electrically  or  mechanically,  or  by  a combination  of  both.  Band- 
width of  the  loops  need  only  be  sufficient  to  ensure  that  positive 
feedback  at  the  dominant  structural  modes  does  not  occur.  Sensitivity 
to  structural  vibrations  within  the  passband  is  reduced  by  notch  filters 
at  the  inputs  to  the  stability-augmentation  system. 

Tentative  characteristics  of  analog  accelerometer  and  rate-gyro 
loops  are  given  in  Tables  2-1  and  2-2.  Both  loops  are  essentially 
second-order  low-pass  filters,  i.e.,  the  transfer  functions  are  of  the 
form 

T = -2 2 ‘2-1) 

s + 2 Cm  s + ui 
n n 


Table  2-1.  Characteristics  of  flight-control 
accelerometers  for  MIRA. 


Parameter 

Value 

Range 

*10  g 

Threshold 

0.01  g 

Natural  Frequency 

15  to  30  Hz 

Damping  Ratio 

0.6  to  0.9 

Linearity 

0.5%  of  full  scale 

Output* 

0.5  V/g 

* Norma]  accelerometer. 
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Table  2-2.  Characteristics  of  flight-control  rate 
gyros  for  MIRA. 


Parameter 

Value 

Range 

±400  deg/s 

Threshold 

±0.1  deg/s 

Natural  Frequency 

15  to  30  Hz 

Damping  Ratio 

0.6 

Linearity 

0.25%  of  full  scale 

Output* 

0.125  V/deg/s 

it 

Roll  rate  gyro. 


The  resulting  time  delay  is  of  interest  when  one  considers  mech- 
anizing a digital  stability-augmentation  or  flight-control  system.  The 
phase  angle,  $(w),  of  the  output  is  given  by 


0 (lu) 


(2-2) 


And  the  time  delay,  of  the  filter  as  a function  of  frequency,  t^Crn), 
is  the  slope  of  the  phase  curve 


td(w) 


dij  (ui) 
do) 


(2-3) 


For  median  values  of  the  filter  parameters,  i.e.,  5 = 0.707  and 

= 130  radians/second,  the  time-delay  curve  is  shown  in  Figure  2-1. 
Over  the  range  of  frequencies  corresponding  to  the  fundamental  struc- 
tural resonances  (6  to  14  hertz),  the  time  delay  is  seen  to  be  12 
to  13  milliseconds. 


Bearing  in  mind  that  the  outputs  of  the  INS  sensors  are  incremen- 
tal velocities  and  angles,  and  hence  need  to  be  differentiated,  the 
above  considerations  would  require  a minimum  sampling  frequency  of 
40  to  50  times  per  second  for  equivalent  performance. 
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Angular  Frequency  (rad/s) 

0 26  52  78  104  130  256  182  208  234  260 


Figure  2-1.  Time  delay  for  second-order  low-pass  filter. 

2.5.4  INS  Requirements  for  Tactical  Applications 

A brief  survey  of  the  INS-related  requirements  for  advanced  tac- 
tical aircraft  displays  a rather  narrow  range  of  acceptable  position 
and  velocity  errors,  but  a much  wider  range  of  tolerable  attitude  and 
attitude  (or  body)  rate  errors.  In  some  cases,  the  accuracy  require- 
ments mandate  the  use  of  position  and  velocity  resets  just  prior  to 
target  acquisition.  In  other  cases,  the  pointing  and  tracking  require- 
ments could  not  be  met  even  by  an  aided  INS,  Disregarding  the  most 
stringent  requirements,  the  ranges  of  acceptable  INS  accuracies  are 
tabulated  in  Table  2-3.  From  the  data  of  Table  2-3,  a set  of  typical 
INS  performance  requirements  is  summarized  in  Table  2-4.  The  indicated 
position-error  rate  corresponds  to  uncompensated  gyro-drift  rates  on 
the  order  of  0.01  degree/hour.  The  latter  would  produce  an  initial 
azimuth  misalignment  of  about  1.0  milliradian,  at  40  degrees  latitude 
with  perfect  accelerometers.  In  a gimbaled  platform,  the  gyro-drift, 
accelerometer,  and  initial-alignment  errors,  which  are  correlated  at 
entry  into  the  navigation  mode,  will  remain  correlated  during  the 
mission.  For  a strapdown  system  this  is  not  the  case,  and  the  error 
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propagation  will  be  strongly  mission-dependent.  Generally,  the  strap- 
down  INS  navigation  and  attitude  errors  will  be  considerably  larger 
than  those  in  a gimbaled  INS  with  the  same  inertial  components  due  to 
the  vastly  greater  angular  rates  encountered.  Gyro  misalignment  and 
scale-factor  errors  are  also  greatly  amplified  by  these  higher  angular 
rates.  In  a strapdown  system,  an  accelerometer  bias  of  100  pg  could 
quite  easily  produce  a peak-velocity  error  in  excess  of  3.0  knots. 

Table  2-3.  INS  accuracy  ranges  for  advanced  tactical  applications. 


Parameter 

Value 

Position 

Horizontal 

0.5  to  1.0  nmi/h  (CEP) 

Vertical ) 

167  to  200  ft 

Velocity ' “ ^ 

Each  Axis 

3.0  to  5.0  kn 

Attitude ^ 

Level 

0.25  to  6.0  mrad 

Heading 

0.8  6 to  20.6  itiin 

Attitude  Rate 

Each  Axis 

0.15  to  1.0  mrad/s 

Body  Rate 

Each  Axis 

0.0086  to  0.057  deg/s 

(1)  May  require  position/velocity  updating  near  target. 

(2)  Inertial  plus  air-data  errors. 

(3)  Not  rms,  but  instantaneous  errors.  For  some  applica- 
tions 1 knot  is  not  enough. 

(4)  Heading  error  is  greater  than  level  errors. 

Table  2-4.  Typical  INS  requirements  from  target-acquisition, 
weapon-delivery,  and  fire-control  requirements. 


Parameter 

Value 

Position  Horizontal 

0.7  nmi/h  (CEP) 

Vertical 

180  ft* 

Velocity  Each  Axis 

4.0  kn** 

Attitude  Level 

4 . 0 rniri** 

Heading 

6 . 0 rniri** 

Attitude  Rate  Each  Axis 

0.02  deg/s** 

* 

Includes  effects  of  nonstandard  lapse  rate. 

* 

lo  values  at  target  acquisition. 


10 


with  these  considerations  in  mind,  existing  and  proposed  strapdown 
systems  and/or  components  were  surveyed  in  an  effort  to  find  systems 
or  components  with  demonstrated  capability  of  meeting  both  the  per- 
formance and  the  dynamic  requirements  of  advanced  tactical  missions. 
Failing  this,  systems  or  components  with  suitable  demonstrated  labora- 
tory performance  in  developmental  quantities  might  be  considered. 

Before  attempting  to  summarize  strapdown-system  performance,  it 
is  useful  to  tabulate  the  more  physical  characteristics  (e.g.,  size, 
weight,  power)  and  cost,  MTBF,  etc.  of  a typical  strapdown  INS  from 
the  survey.  These  results  are  given  in  Table  2-5. 

Table  2-5.  Typical  strapdown  INS  characteristics. 


Parameter 

Value 

Volume 

<1000  in.^ 

Weight 

<30  lb 

Power  - Peak 

<250  W 

Operating 

<100  W 

Cost  in  Production 
Quantities 

<$30,000 

MTBF 

>3000  h 

Acceleration 

>10  g 

Angular  Rate 

>300  deg/s 

From  the  characteristics  in  Table  2-5,  it  is  seen  that  the  physi- 
cal and  cost  characteristics  of  the  typical  strapdown  INS  are  generally 
compatible  with  the  intended  application.  However,  the  level  of  per- 
formance is  very  environment  sensitive,  and  ranges  from  0.6  to  4.0 
nautical  miles/hour  (CEP).  Furthermore,  the  initial-alignment  method 
and  reaction  time  are  frequently  unspecified.  It  should  be  noted  that 
the  reaction  time  for  tactical  applications  is  normally  10  minutes  or 
less . 

Table  2-5  is  representative  of  the  future  class  of  strapdown  sys- 
tems comprising  floated  rate  integrating  gyros,  tuned  rotor  gyros,  laser 
gyros,  ESGs , and  even  multisensors.  It  is  representative  of  proposed 
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systems  and  those  of  which  a limited  number  have  been  built.  Each  of 
these  sensors  has  its  claimed  advantages  and  disadvantages,  but  it  is 
clear  that  the  multisensors  would  offer  the  minimum  number  of  components 
required  to  perform  three-axis  rate  and  specif ic- force  measurement. 

For  the  purposes  of  demonstrating  the  MIRA  simulation,  however,  laser 
gyros  were  chosen  as  the  initial  rate  sensors  to  be  modeled,  and 
single-degree-of-freedom  accelerometers  as  the  specific-force  measuring 
devices . 

The  incremental  velocity/linear  acceleration  sensors,  in  the 
strapdown  systems  surveyed,  are  generally  of  the  pendulous  single-axis 
type,  with  an  analog  rebalance  loop  and  subsequent  conversion  to  digital 
format  (either  direct  or  voltage-to-pulse  rate  and  summation) . The 
most  common  accelerometer  is  the  Sunstrand  (Q-Flex) , which  is  a quartz 
flexure  hinged-pendulum  device  with  adequate  range  and  performance. 

One  disadvantage  of  trying  to  recover  instantaneous  acceleration  from 
the  incremental  velocity  (digital  output)  of  the  accelerometer  is  the 
fact  that  the  quantization  tends  to  be  on  the  order  of  0.032  foot/ 
second/pulse  or  1,000  pulses/second/g  which  militates  against  the  re- 
covery of  instantaneous  acceleration,  i.e.,  with  acceptably  small  time 
delay.  However,  separate  processing  of  the  analog  output  could  yield 
the  desired  acceleration  with  negligible  delay. 

2.2.5  Minimum  Sensor  Configuration 

It  is  expected  that  during  the  1980  time  frame,  triply  redundant 
sensor  systems  will  be  standard  on  advanced  tactical  fighters.  There- 
fore, the  minimum  sensor  configuration  would  consist  of  three  strapdown 
systems,  each  located  at  different  fuselage  stations  in  the  aircraft. 
These  would  be  FS  77,  FS  186,  and  FS  383  for  the  simulation  vehicle, 
which  are  considered  the  optimum  locations  for  the  normal  accelerom- 
eter and  roll  rate  gyro,  the  lateral  accelerometer,  and  the  pitch  and 
yaw  rate  gyros  respectively.  Thus,  when  all  components  are  operating, 
there  would  be  one  stability-augmentation  sensor  of  each  kind  in  its 
optimum  location  (master)  plus  two  backup  sensors  in  nonoptimum  loca- 
tions. Presumably,  filtering  would  be  applied  to  the  outputs  of  the 
sensors;  at  least  those  in  the  nonoptimum  locations  to  minimize  the 
differences  in  indicated  outputs  relative  to  the  masters.  Lever-arm 
corrections  would  certainly  have  to  be  applied  for  at  least  two 
thirds  of  the  sensors. 
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At  each  of  the  three  locations,  a complete  strapdown  IMU  (and 
possibly  its  own  navigation  computer)  would  be  mounted.  Preprocessing 
of  the  stability-augmentation  signals  would  also  be  performed  at  these 
locations  either  in  an  existing  navigation  computer  or  separate  micro- 
processor (s ) . 

The  signals  would  be  mixed  and  the  stability-augmentation/flight- 
control  computations  would  be  performed  in  one  (or  more)  centralized 
computer (s),  which  would  provide  the  actuator  signals. 

Radio  navigation  aids,  air  data  (with  the  possible  exception  of 
barometric  altitude) , and  the  navigation  and  attitude  data  from  the 
three  strapdown  INS  would  be  processed  in  the  same  (or  other)  computer (s) 
to  provide  "best”  estimates  of  the  navigation-data  and  navigation- 
parameter  updates.  The  various  pointing,  tracking,  fire-control, 
weapon-delivery,  and  display  computations  would  likely  be  performed 
in  this  computer  complex  also. 
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SECTION  3 


MIRA  SIMULATION  DEVELOPMENT 


3 . 1 Introd action 

A major  goal  of  the  present  program  is  the  development  of  a simu- 
lation capability  to  evaluate  the  MIRA  concept.  In  addition,  the  capa- 
bility of  the  simulation  is  to  be  demonstrated.  This  section  describes 
the  simulation,  and  details  the  assumptions  made  in  generating  it,  as 
well  as  the  models,  equations,  and  data  programmed.  A presentation  of 
the  data  obtained  from  the  use  of  the  simulator  and  a discussion  of  the 
corresponding  results  follows  in  Section  4. 

The  simulation  goals  of  the  MIRA  project  can  be  achieved  in 
several  different  ways.  A series  of  simulator  goals  have  been  es- 
tablished which  will  provide  the  maximum  long-term  benefit  to  the 
program  and  concept  evaluation.  One  of  these  goals  is  the  develop- 
ment of  a real-time  simulation  which  will  lessen  computer  operating 
costs  and  allow  human  interaction  with  the  system.  A fixed-base 
cockpit  has  been  tied  into  the  system  allowing  pilot  evaluation  of 
the  candidate  MIRA  configurations  prior  to  flight  test.  Another  CSDL 
simulation  goal  is  the  development  of  an  all-digital  simulation  capa- 
bility. There  are  many  factors  which  make  an  all-digital  simulation 
most  attractive  when  compared  to  a system  utilizing  analog  computers. 
Among  them  are  ease  of  setup,  repeatability,  reliability,  etc.  A third 
feature,  incorporated  into  the  MIRA  simulation,  is  flexibility.  Pro- 
visions have  been,  and  will  continue  to  be,  included  in  the  simulation 
which  allows  ease  of  change  of  parameters,  sensor  configurations, 
sensor  locations,  program  subroutines  and  models,  aircraft  aerodynam- 
ics, and  (in  general)  almost  any  portion  of  the  simulation.  The  long- 
term benefit  of  each  of  these  goals  to  the  MIRA  program  is  clearly 
evident. 

The  MIRA  simulation  utilizes  a Xerox  Data  System  Model  (XDS)  9300 
general-purpose  digital  computer  as  the  basic  computation  unit.  This 
machine  has  a 32,768-word  magnetic  core  memory,  24-bit  word  length,  and 
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1.75-microsecond  machine  cycle  time.  In  the  near  future,  this  computa- 
tion capability  will  be  expanded  by  interfacing  two  Honeywell  124-series 
digital  computers  with  the  XDS  9300  to  operate  as  parallel  processors. 

One  of  these  units  is  a basic  124  expanded  to  a Honeywell  224  memory 
configuration  of  16k  worus,  while  the  other  unit  is  a basic  124  with 
8k  words  of  memory. 

The  basic  program  plan  calls  for  the  development  of  a non-real- 
time version  of  the  MIRA  simulation,  programmed  in  Fortran  on  the 
XDS  9300,  for  purposes  of  checkout  and  to  obtain  preliminary  data  and 
results.  This  same  simulation  capability  is  to  be  expanded  to  real-time 
operation  using  fixed-point  assembly  language.  The  124  computers  are 
needed  for  this  phase  of  the  operation.  At  present,  the  evolution  of 
the  non-real-time  simulation  has  been  completed.  The  conversion  of 
this  simulation  to  real-time  operation  has  begun,  but  only  includes 
the  rigid-body  equations  of  motion  of  the  vehicle  at  this  point.  The 
next  step  in  the  program  involves  bringing  the  real-time  simulation 
capability  up  to  the  level  of  the  non-real-time  simulation. 

3 . 2 General  Description  of  MIRA  Simulation  Characteristics 

A block  diagram  of  the  system  programmed  in  the  non-real-time 
simulation  is  presented  in  Figure  3-1.  The  nonlinear  six-degree-of- 
freedom  equations  of  motion  of  the  vehicle  have  been  programmed.  A 
spherical  rotating  earth  model  is  used.  The  aerodynamics  incorporated 
into  the  simulation  are  valid  for  a cleanly  configured  vehicle  in 
subsonic  flight.  More  specifically,  the  aerodynamics  are  programmed 
as  nonlinear  polynomial  functions  valid  over  an  envelope  of  +15  degrees 
in  sideslip  angle,  and  up  to  30  degrees  in  angle  of  attack.  This  large 
envelope  allows  a great  deal  of  flexibility  in  the  use  of  the  simulation, 
without  any  extensive  reprogramming. 

Three  longitudinal  and  three  lateral-directional  structural  modes 
have  been  added  to  the  rigid-body  equations  of  motion.  A stability  and 
control  augmentation  system  has  also  been  included  to  improve  the  han- 
dling qualities  of  the  vehicle.  Structural-mode  filters  are  also 
included  in  the  augmentation  system.  It  is  possible  to  enter  time- 
history  profiles  into  the  simulation  to  generate  typical  aircraft  maneu- 
vers. The  effects  of  winds  and  turbulence  on  the  vehicle  are  another 
feature  of  the  simulation. 
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Figure  3-1.  MIRA  simulation  block  diagram. 

Models  of  flight-control  sensors,  and  navigation-system  laser 
rate  gyros  and  quartz-flexure  acceleron eters  are  i ncor])orated  into  the 
MIRA  simulation.  It  is  possible  to  independently  locale  each  of  tl-.ese 
sensors,  and  to  use  the  navigation  sensors  for  flight  control.  An 
onboard  Local  Vertical  Wander  Azimuth  (L\AVA)  navigation-system  model 
has  also  been  programmed.  This  model  was  obtained  from  the  Standard- 
ized Software  Development  Program.  Errors  between  vehicle  and 

navigation-system  parameters  are  also  generated  for  the  purpose  of 
evaluating  different  MIRA  configurations.  A detailed  discussion  of 
each  of  these  aspects  of  the  MIRA  simulation  is  presented  in  the  follow- 
ing sections  of  the  report  so  that  a better  understanding  of  its  overall 
capability  can  be  obtained. 

One  flight  condition  has  been  selected  for  development,  demonstra- 
tion, and  evaluation  of  the  MIRA  simulation  capability.  The  vehicle  is 
at  an  altitude  of  5,000  feet  and  Mach  0.5.  The  trim  velocity  is  549 
feet/second,  and  the  trim  angle  of  attack  is  2.06  degrees.  Selection 
of  this  point  in  the  flight  envelope  was  dictated  primarily  by  the  fact 
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that  consistent  data  are  available  for  each  portion  of  the  simulation 
at  this  flight  condition,  while  simultaneously  allowing  the  goals  of 
the  MIRA  program  to  be  satisfied.  The  availability  of  time  histories, 
vehicle-model  parameters,  and  other  data  useful  for  system  checkout, 
also  dictated  the  choice  of  this  flight  condition  as  the  nominal  one. 

The  reader  is  cautioned  that  some  of  the  data  given  in  this  report 
apply  specifically  to  the  nominal  flight  condition,  and  will  have  to 
be  modified  should  another  flight  condition  be  considered. 

3 . 3 Vehicle  Six-Degree-of-Freedom  Nonlinear  Equations  of  Motion 

This  section  contains  a development  of  the  six-degree-of-f reedom 
nonlinear  equations  of  motion  of  the  vehicle.  The  development  closely 
follows  that  given  in  Reference  6. 

Four  orthogonal-axes  systems  must  be  defined  to  develop  the  appro- 
priate equations  of  vehicle  motion.  They  are  depicted  in  Figures  3-2, 
3-3,  and  3-4,  and  are  defined  as: 

(1)  The  inertial  frame  which  is  taken  as  fixed  in  space,  and 
for  which  Newton's  Laws  of  Motion  are  valid. 

(2)  An  earth-centered  frame  which  rotates  with  the  earth. 

(3)  An  earth-surface  frame  which  is  parallel  to  the  earth's 
surface,  and  whose  origin  is  at  the  vehicle's  center  of 
gravity  (eg)  defined  in  north,  east,  and  down  directions. 
This  frame  moves  with  the  vehicle,  and  its  center  is 
defined  with  respect  to  the  earth-centered  frame  by  the 
longitude,  latitude,  and  altitude  of  the  vehicle. 

(4)  The  conventional  body  axes  are  selected  to  represent  the 
vehicle.  The  center  of  this  frame  is  at  the  eg  of  the 
vehicle  and  its  components  are  forward,  out  of  the  right 
wing,  and  down. 

Several  major  assumptions  have  been  made  in  the  basic  development 
of  the  vehicle  equations  of  motion.  These  assumptions  are; 

(1)  The  vehicle  equations  are  written  in  the  body-axes  coordi- 
nate frame. 

(2)  A spherical  earth  rotating  at  a constant  angular  velocity 
is  assumed. 

(3)  The  vehicle  aerodynamics  are  nonlinear. 
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Figure  3-4.  Orientation  angles. 


(4)  The  iindisturbed  atmosphere  rotates  with  the  earth. 

(5)  Tile  winds  arc  defined  v/ith  respect  to  the  earth. 

(6)  The  gusts  are  defined  with  respect  to  the  vehicle. 

(7)  An  inverse-sc;uare  gravitational  law  is  used  for  the 
spherical-earth  model. 

(8)  The  mass  of  the  vehicle  is  constant. 

(9)  Flight  over  the  poles  is  excluded  due  to  the  use  of  lati- 
tude as  a variable. 

(10)  The  gradients  of  the  lovvr-frequency  winds  are  small  enough 
to  be  neglected. 

Other  assumptions  which  have  been  made  during  the  development  of  the 
MIRA  simulation  are  mentioned  as  they  arise  during  the  discussion. 

The  force  equations  of  the  vehicle,  with  respect  to  the  inertial 
reference,  are  given  in  vector  form  as: 

IVI  = i [FI  - [ , ) • |Vj  (3-1) 

while  the  moment  equations  are  written  as 

III  h J [Ml  - [wxl  [I]  |..l  (3-2) 
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where 
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Expansion  of  these  vector  equations  leads  to  the  following  familiar 
scalar  equations. 


u = — + rv  - qw 

m 


V = 

m 


ru  + pw 


w = — + qu  - pv 


P = f ["xz^  ^ - ^z>  ^ P^^xz  ^ \]  '3-7) 

^ ^ [^xz'*^'  - P'>  P>^<3z  - 3x>  " \]  '3-8) 

^ f ["xzP  - ^ '3x  - 3y)Pq  + M,]  (3-9) 

z 

The  vehicle  force  vector  may  be  further  defined  in  terms  of  the 
aerodynamic  structural  mode  and  gravitational  forces,  and  the  thrust  of 
the  right  and  left  engines  as 


(F)  = [F 


aero^  '^grav^  ^ 


(3-10) 


The  aerodynamic  force  vector  is  defined  as 


[F  ] = qS  C 

aero^  ^ y 


(3-11) 
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where  *3  ~ J gravitational  force  vector  is  given  by 
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An  inverse-square  law  is  used  to  compute  g,  i.e., 


(3-13) 


The  force  components  due  to  the  engine  thrust  and  the  structural  modes 
are  developed  in  more  detail  in  a latter  section.  Expansion  of  the 
force  vector  results  in  the  following  scalar  equations. 
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The  vehicle  moment  vector  may  also  be  defined  in  terms  of  the 
aerodynamic  structural  mode  and  thrust  moments  as 


[M] 
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[M  ] = qS 
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where 


Acg  = eg  - 0.33 
^ ^pos 


(3-19) 


The  incremental  terms,  C Acg  subtracted  from  C and  C Acgr-  added  to  C , 
account  for  the  location  of  the  aircraft  eg  at  a position  different 
from  the  nominal  value  of  0.33c,  at  which  the  aerodynamic  data  are 
defined.  A development  of  the  thrust  and  structural-mode  moments  is 
detailed  in  a later  section  of  this  discussion.  However,  the  moment 
vector  can  be  expanded  in  terms  of  tlie  components  of  these  vectors  as 
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(3-21) 


My  = qsc(c^  - C^Acg)  + 
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The  aerodynamic  forces  and  moments  are  computed  using  the  total 
velocity  of  the  aircraft  eg  relative  to  the  local  atmosphere,  which  is 
given  by 


V, 

rt 


(3-23) 


Tlie  vector  equation  for  the  body-axes  velocity  with  respect  to  the 
local  atmosphere  is 
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In  terms  of  the  x,  y,  and  z components  of  each  of  these  vectors,  the 
atmosphere  relative-velocity  components  are  defined  by 
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where  ' '^rel'  '^rel  linear  body-axes  components  of  the 
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The  winds  are  defined  with  respect  to  the  earth  with  components 
emanating  from  the  north,  east,  and  below  the  vehicle.  Thus,  a north 
wind  blows  from  north  to  south,  etc.  The  body-axes  wind  com;^onents  are 
tiien  given  by 
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Angle  of  attack  and  sideslip  angle  are  used  in  the  definition  of 
the  vehicle  aerodynamics.  These  variables  are  given  by  the  expressions 


B = 


(3-33) 


(3-34) 


T'~a  angular  velocities  of  the  vehicle  with  respect  to  the  local 
atmosphere  also  have  to  be  computed  for  tlie  vehicle  aerodynamics.  They, 
however,  do  not  turn  out  to  be  as  complicated  as  the  linear  velocities 
with  respect  to  the  local  atmosphere.  The  assumption  that  the  spatial 
gradients  of  the  low-frequency  winds  are  small  enough  to  be  neglected, 
implies  that  the  angular  components  of  the  wind  can  be  neglected.  A 
similar  assumption  is  made  regarding  the  angular  velocities  of  the 
gusts.  The  overall  result  is  that  the  angular  velocities  of  the  body 
axes  with  respect  to  the  atmosphere  are  the  angular  velocities  of  the 
body  axes  with  respect  to  the  earth.  They  are  given  by  the  angular 
velocities  of  the  body  axes  with  respect  to  inertial  space,  less  the 
angular  velocities  of  the  earth-centered  frame  with  respect  to  the 
inertial  reference  in  the  body-axes  coordinate  system.  In  equation 
form,  the  result  is 
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It  is  necessary  to  compute  the  variable  a,  since  the  coefficient 
is  present  in  the  aerodynamics  of  the  vehicle,  and  strongly  affects 
its  short-period  response.  Several  techniques  can  be  used  to  obtain  a 
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One  is  by  back-differencing  a.  This  is  an  unacceptable  method  when 
high-frequency  components  are  present  in  the  a response,  such  as  due 
to  structural  modes  and/or  gusts.  An  alternate  technique  is  to  compute 
a from  the  expression 


(3-38) 


''''A  - 

''I 

Z X 


(3-39) 


Use  of  tliis  latter  expression  presents  an  additional  i,roblem  in  that  it 

necessitates  the  computation  cf  and  . These  latter  variables 

^x  z 

can  be  obtained  from  the  vector  differentiation  of  Eq . (3-24).  However, 

fairly  complicated  expressions  result.  Simplifications  to  the  true 
expressions  for  and  are  discussed  in  Reference  6.  It  is  shoivn 

that  an  acce;ptably  accurate  solution  for  an.J  may  1;;:  obtained 

X z 

if  it  is  assumed  that  the  earth  is  rotating  slowly,  and  that  these 
acceleration  variables  are  used  only  to  compute  a and  not  integrated 
to  obtain  , or  a.  Neither  of  these  assumptions  is  prohibitive 

in  this  particular  case. 

The  expression  resulting  from  the  vector  differentiation  of 
Eq.  (3-24),  and  the  incorporation  of  the  assumptions  just  discussed  is 

[V^l  = [V]  + [AelXltTEsel'^  {[V^]  + [W]  - (3-40) 

It  is  indicated  in  Reference  6 that  the  effect  of  the  gusts  on  the 
terms  of  the  previous  expression  is  only  on  the  order  of  10  percent, 
and  that  neglecting  will  not  cause  serious  errors.  This  assump 

tion  is  made  in  the  MIRA  simulation  since  it  is  desired  to  avoid  the 
differentiation  of  the  gust  signals.  The  matrix  equation  for  the  co;n- 
ponents  of  [^^1  are  then  given  by 
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The  specific  expressions  programmed  are 

• • 

^ = u 

X 

• • 

V,  = w 


“ei 

■ u . , "I 

wind 

'^EI 

+ 

'^wind 

'^EI- 

- '^wind  - 

(3-42) 
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^'^EI  ■*'  %ind' 

+ Prel  <''ei 

^ ''wind' 

(3-44) 

The  motions  of  the  vehicle  relative  to  the  earth  (the  trajectory 
kinematics)  are  described  by  the  motion  of  the  earth-surface  axes.  The 
earth  relative  velocities  are  given  by 


"x ' 
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r , -1 
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^'^esb' 
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(3-45) 
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Making  use  of  Eq.  (3-30)  leads  to  an  alternate  definition  in  terms  of 
the  vehicle  velocities  relative  to  earth,  which  is 


-1 


w. 


rel 
■'rel 
rel  J 


(3-47) 


The  position  of  the  vehicle  is  described  by  its  latitude,  longitude, 
and  altitude.  These  quantities  are  obtained  by  integrating  the  earth 
velocities  of  the  vehicle  in  spherical  coordinates.  The  appropriate 
equations  for  integration  are 


h = -Z 


^ ^ 


(3-48) 

(3-49) 


o = 


Rg  cos  A 


An  auxiliary  variable,  the  flightpath  angle  y,  is  defined  as 

h 


7 = sin“^ 


rel 


(3-50) 


(3-51) 
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Euler  angles  defining  the  earth  surface  to  body-axes  trans- 
are  defined  in  Figure  3-4  and  obtained  by  integration  of  the 
equations . 


0 = 

10^  + 0 sin  0 

(3-53) 
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= U' 

y 

cos  0 - oj^  sin 
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(3-54) 
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(3-55) 

cos  9 z cos 

e 

The  angular  velocity  components  of  the  vehicle  about  its  body  axes  with 
respect  to  the  earth-surface  axes  system  co^ , and  u>^)  are  given  by 


rto  1 
X 

■Prel' 

a cos  A 

LO 
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^rel 

^'^esb' 

-A 
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L z J 

. ^rel. 

-a  sin  A 

This  expression  is  nothing  more  than  the  angular  body  rates  of  the 
vehicle  with  respect  to  the  earth,  less  the  angular  rates  of  the  earth- 
surface  axes  system  with  respect  to  the  earth  in  body  axes. 


The  transformatiori  from  the  earth-surface  axes  system  to  the 
body-axes  frame  is  achieved  by  a yaw,  pitch,  and  roll  about  the  longi- 
tudinal, lateral,  and  normal  axes,  respectively.  The  resultant  matrix 
is 


cos  6 cos  Ip 

sin  ; sin  6 cos  i)/ 

- cos  (?  sin  i|) 

cos  ((>  sin  6 cos  ip 

+ sin  ((>  sin  ip 


cos  9 sin 
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+ cos  0 cos  4 

cos  0 sin  6 sin  ip 

- sin  0 cos  0 


-sin  0 
sin  0 cos  e 

cos  0 cos  0 

(3-57) 


3 . 4 Vehicle  Aerodynamics 

Expressions  for  the  aerodynamic  coefficients  of  the  vehicle  are 
presented  in  Reference  1 in  polynomial  form.  These  polynomials  are 
functions  of  a,  a,  6,  V^,  the  body-axes  angular  rates  relative  to  the 
atmosphere,  and  the  aerodynamic  control  surfaces.  The  control  sur- 
faces consist  of  a combined  spoiler  and  aileron  deflection  with 

1.5  degrees  of  spoiler  deflection  for  each  degree  of  aileron  deflection 
(6g) , the  horizontal  tail  or  stabilator  (6^) , and  the  rudder  (6^) . 
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Expressions  are  available  for  three  angle-of-attack  ranges  (a  £ 15  de- 
grees, 15  degrees  a £ 30  degrees,  and  30  degrees  a < 45  degrees), 
since  the  intent  of  the  study was  to  identify  the  high  angle-of- 
attack  aerodynamic  coefficients  of  the  vehicle.  The  first  two  sets 
of  these  polynomials  have  been  programmed  into  the  MIRA  simulation  since 
the  vehicle  stalls  at  an  angle  of  attack  of  approximately  23  degrees, 
and  a definition  beyond  this  boundary  is  felt  to  be  unnecessary  at  this 
time.  The  aerodynamic  expressions  follow. 


i-range:  a s 15'^ 
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C = -6.61  X 10  ^ - 2.67  X lo“^a  - 6.48  x 10 
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Tliese  data  are  valid  only  for  a vehicle  in  the  clean  configuration;  land- 
ing gear,  flap  effects,  etc.,  are  not  accounted  for.  In  addition,  the 
aerodynamic  model  ife  valid  only  for  subsonic  flight  with  M = 0.7  about 
the  upper  limit. 

A few  words  of  caution  are  in  order  concerning  the  use  of  these 
polynomials  defining  the  aerodynamics  of  the  vehicle.  The  transition 
between  the  angle-of-attack  ranges  is  not  smooth,  and  discontinuities 
are  present  at  the  boundaries  for  some  of  the  variables.  The  time 
histories  obtained  for  those  cases  v^here  a transition  from  one  angle-of- 
attac)c  range  to  another  exists  may  exhibit  discontinuities.  The  impli- 
cation of  these  unwanted  and  unrealistic  responses  should  be  investi- 
gated to  assess  their  effect  on  the  validity  of  the  study  being 
conducted. 
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. 5 Thrust  Forces  and  Moments 


Each  of  the  force  and  moment  equations  contains  thrust  terms  due 
to  the  fact  that  the  engines  are  located  on  the  sides  of  the  fuselage 
and  are  inclined  at  an  angle  with  respect  to  the  vehicle  body  axes. 

As  discussed  in  Reference  1,  the  orientation  of  the  thrust  is  defined 
by  two  angles,  a,j,  and  By  definition,  is  the  angle  between  the 

thrust  vector  and  a line  parallel  to  the  x body  axis,  measured  in  a 
plane  parallel  to  the  x axis  and  containing  the  thrust  vector;  and 
is  the  angle  between  the  x-y  plane  and  the  plane  parallel  to  the  x axis, 
and  containing  the  thrust  vector.  For  the  selected  model,  X^  = 5.25  de- 
grees, and  = 87.28  degrees.  Solving  for  the  thrust  forces  along  the 
body  axes  by  resolving  vector  components  leads  to  the  following  expressions. 


Tj^  cos  X^ 
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The  thrust  moments  are  generated  because  the  point  of  thrust  appli- 
cation is  offset  from  the  vehicle  eg.  For  the  right  engine,  the  point 
of  application  is  defined  by  the  following  body-axes  moment  arms. 


.608  + 

(3-72) 

> 

II 

1.985  ft 

(3-73) 

-0.391  ft 

(3-74) 

That  is,  the  point  of  application  is  0.608  foot  in  front  of  0.33  c along 
the  X body  axis,  1.985  feet  out  the  right  wing  along  the  y body  axis, 
and  0.391  foot  above  the  x-y  plane  along  the  z body  axis.  The  increment 
(cg^^^  - 0.33)c  in  the  Ax^  term  accounts  for  eg  location.  For  the  left 

engine,  the  appropriate  moment  arms  are 


.608  + 

(eg 

^pos 

0.33)c  ft 

(3-75) 
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The  expressions  for  the  thrust  moments  using  these  moment  arms,  and  the 
expressions  for  the  thrust  forces  given  previously  are 
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3 . 6 Trim  Algorithm 

The  MIRA  simulation  includes  an  automatic  trim  routine  in  which 
a,  6,  0,  and  the  control  deflections  are  calculated  so  that  the 
vehicle  body-axes  accelerations  and  side  and  normal  velocities  with 
respect  to  the  earth  surface  are  ideally  zero  or,  in  the  practical 
case,  very  small.  An  iterative  scheme  is  used  to  adjust  these  param- 
eters to  the  required  trim  values.  Initialization  of  each  of  the  sub- 
routines of  the  simulation  also  occurs  during  the  automatic  trim  mode. 
The  vehicle  is  always  trimmed  with  wings  level,  so  that  if  a crosswind 
is  present  the  nose  of  the  vehicle  is  turned  into  the  wind  at  trim. 

The  section  includes  a brief  discussion  of  the  trim  algorithm 
included  in  the  MIRA  simulation  to  present  the  basic  ideas  which  it 
incorporates.  A discussion  of  the  initialization  aspect  of  the  simula- 
tion is  not  included  since  this  is  felt  to  be  unnecessary,  in  addition 
to  being  very  unwieldy. 

A desired  altitude  and  Mach  number  are  selected  for  the  vehicle, 

as  well  as  initial  estimates  for  the  control  deflections,  a and  6. 

Then,  V is  computed  from  the  speed  of  sound  and  Mach  number  by 
^t 

V-  = V • M (3-80) 

A^  s 

, and  are  then  calculated. 
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V,  tan  a. 
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■t  ( ''itj 


(3-81) 
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Tlie  earth-relative  linear  velocities  are  then  calculated  from  the  afore- 
mentioned parameters  using  the  equations 


rel. 


rel. 


rel. 
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V + w • j w 
A wind  q 

^t 


(3-84) 
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and  ‘^rel'  ^rel  fixed  at  zero. 

The  trim  pitch  angle  is  chosen  so  that  Z = 0.0  by  the  formula 


9^  = tan 


"fe) 


(3-87) 


and  the  trim  heading  angle  is  chosen  so  that  Y = 0.0  by  the  formula 

-1 


= -tan 


fe) 


(3-88) 


Both  of  these  expressions  assume  wings  level  flight  so  that  = 0.0. 
Furthermore,  the  expression  for  assumes  that  cos  9^  = 1. 

The  vehicle  equations  of  motion  are  used  to  calculate  w,  and 
is  adjusted  so  that  w^  = 0.0.  An  iterative  scheme  is  used 

+ 0.0003W 


a = a 

k k-1 


(3-89) 


until  w^  is  sufficiently  small.  A similar  scheme  is  used  to  zero  v^ 
by  adjusting  6^  according  to  the  relationship 
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, V , and  V are  then  recalculated  using  Eq . (3-81),  (3-82), 

^t  ^t  ^t 

and  (3-83) . 

The  control  deflections  are  now  adjusted  to  null  the  other  body- 
axis  accelerations,  and  u^  is  zeroed  by  adjustment  of  each  throttle 
using  the  expression 


k-1 

O.lil 

(3-91) 
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With  the  completion  of  this  phase  of  the  algorithm,  control  reverts 
back  to  the  beginning  of  the  trim  routine,  and  the  whole  cycle  continu- 
ously repeats  until  the  initiation  of  the  compute  mode. 


3 . 7 Structural  Modes 

The  equations  of  motion  include  the  effects  of  three  longi- 
tudinal and  three  lateral  directional  aeroelastic  modes.  References  2 
and  3 contain  the  appropriate  equations  and  coefficients  valid  for 
small  perturbations  of  the  aircraft  from  straight  and  level  flight. 

Consider  the  longitudinal  modes  first.  Modes  involving  strong 
fuselage-stabilator  coupling  were  selected  because  they  are  excited  by 
the  stabilator-control  inputs,  and  contribute  significantly  to  the 
sensor  responses.  Additionally,  vibration  modes  were  selected  which 
adequately  define  the  stabilatorstatic  (as  well  as  dynamic)  aerodynamic 
stability  characteristics.  The  three  longitudinal  modes  accounted  for 
are ; 

(1)  A stabilator-bending  mode,  denoted  by  the  variable  which 

has  a nominal  frequency  in  the  vicinity  of  71.0  radians/ 
second,  and  a nominal  damping  ratio  of  0.012. 

(2)  The  first  vertical-bending  mode,  designated  by  the  variable 

with  a nominal  frequency  of  approximately  87.5  radians/ 
second  and  a nominal  damping  ratio  of  about  0.021. 
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(3)  A stabilator-rotation  mode,  denoted  by  the  variable  at 
a nominal  frequency  of  about  141.0  radians/second  and  a 
nominal  damping  ratio  of  0.058, 

Each  of  these  modes  is  defined  by  a second-order  differential  equation 
with  additional  terms  which  couple  in  the  basic  airframe  responses,  the 
other  modes,  and  the  stabilator  deflections.  In  addition,  the  effect 
of  the  modes  on  the  pitching-moment  and  lift  equations  is  available. 

The  effect  of  the  structural  modes  on  the  longitudinal  equation  of 
motion  is  not  considered  in  References  2 and  3 and,  therefore,  is  not 
considered  in  the  present  simulation  mechanization. 

The  longitudinal  structural-mode  data  are  available  in  terms  of 
dimensional  derivatives.  For  the  MIRA  simulation,  the  coefficients  of 
the  aeroelastic  equations  have  been  nondimensionalized,  and  the  equa- 
tions programmed  with  constant  coefficients  as  a function  of  dynamic 
pressure  and  true  airspeed.  In  addition,  1^  has  been  factored  out  of 
the  aeroelastic  terms  of  the  pitching-moment  equation,  and  m has  been 
factored  out  of  the  aeroelastic  terms  of  the  lift  equation.  This  has 
been  done  for  two  reasons.  First,  the  coefficients  of  the  structural- 
mode terms  in  the  pitching-moment  and  lift  equations  are  strongly 
affected  by  changes  in  1^  and  m,  respectively.  Nondimensionalizing  the 
coefficients  provides  a means  of  compensating  for  changes  in  these 
parameters  easily  and  properly.  Second,  the  simulation  is  more  flexible 
and  allows  the  investigation  of  a wider  flight  regime.  While  it  is 
realized  that  some  of  the  nondimensional  coefficients  of  the  structural- 
mode equations  change  with  flight  condition,  they  are  also  not  available 
for  all  flight  conditions.  It  may  be  possible  to  obtain  an  indication 
of  the  true  results  in  these  latter  cases  due  to  the  nondimensionaliza- 
tion. 

The  longitudinal  rigid-body  modes  and  the  fuselage  and  stabila- 
tor elastic  modes  are  widely  separated  in  frequency.  For  this  reason, 
the  static  effects  of  the  elastic  modes  on  the  rigid-body  coefficients 
are  negligible,  and  the  latter  parameters  can  be  used  without  modi- 
fication. Furthermore,  changes  in  the  vehicle  eg  do  not  affect  the 
structural-mode  coefficients  significantly;  it  is  sufficient  to  ac- 
count for  eg  variations  in  the  rigid-body  coefficients  as  indicated  in 
Section  3.3. 

The  equations  given  in  References  2 and  3 are  in  the  stability- 
axes  coordinate  system.  At  the  nominal  flight  condition  selected  for 
the  MIRA  simulation,  the  pertinent  coefficients  of  these  equations 
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have  been  used  directly  in  the  body-axes  equation  of  the  simulator. 
Since  the  trim  angle  of  attack  at  this  flight  condition  is  about 
2 degrees,  the  effect  is  insignificant. 

The  resultant  longitudinal  structural-mode  equations  programmed 
in  the  MIRA  simulation  are: 


n\  = -1.33ni  - 4436. 46n,  - 22.806* 

J.  i I e 

- q jl.gSrij^  + 2.56ri2  + 13.60ri3  + 239.43A6g| 

- ^ |l0.90AV^  + 321.22qj.g3  + O.eOOn^  + 0.500^2  + 

+ 50.905g  I (3-95) 
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z I 2 e 
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It  is  noted  that  the  aeroelastic  equations  are  a function  of 

• ••  t 

6 and  6 . In  the  first  mechanization  of  these  equations,  6^  was 
e e e 

obtained  from  a first-order  lag  present  in  the  flight-control  system. 
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and  6 obtained  by  back-differencing  .5  . It  was  determined  that  the 

• • • ^ 

amplitudes  of  and  6^  achieved  in  this  fashion  were  inversely  pro- 
portional to  the  sampling  interval,  and  very  sharply  peaked  for  step 
inputs.  Unrealistic  results  were  obtained.  A more  realistic  simulation 

has  been  achieved  by  lagging  6 , and  then  back-differencing  the  lagged 

• • • ® 

to  obtain  6^.  The  time  constant  of  the  lag  is  50  milliseconds.  The 
responses  are  smoother  and  much  less  dependent  upon  the  computation 
interval . 

The  lateral-directional  aeroelastic  modes  are  now  considered. 

The  modes  judged  to  be  most  significant  in  References  2 and  3,  and  for 
which  data  are  available  are: 

(1)  The  fuselage  first  torsion  mode,  designated  by  the  variable 
r,^,  with  a nominal  , frequency  of  39."'  radians/second  and  a 
nominal  damping  ratio  of  0.011. 

(2)  The  wing  first  asymmetric-bending  mode,  designated  by  the 
variable  n^,  with  a nominal  frequency  of  67.6  radians/ 
second  and  a nominal  damping  ratio  of  0.0165. 

(3)  The  fuselage  first  lateral-bending  mode,  designated  by  the 
variable  n,,  with  a nominal  frequency  of  80.8  radians/ 

b 

second  and  a nominal  damping  ratio  of  0.033. 

Each  of  these  modes  is  represented  by  a coupled  linear  second-order 
differential  equation.  However,  unlike  the  longitudinal  modes,  there 
are  no  aeroelastic  terms  in  the  rigid-body  equations. 

In  general,  the  same  approach  and  assumptions  which  applied  to 
the  longitudinal-directional  mode  equations  were  also  applied  to  the 
lateral-directional  mode  equations  (e.g.,  the  equations  are  in  non- 
dimensional  form) . The  three  lateral-directional  modes  do  not  include 
all  of  the  airframe  elastic  modes  through  the  control-system  cutoff 
frequency  needed  to  define  the  flexibility  of  the  primary  surfaces  and 
control  surfaces.  In  References  2 and  3,  this  fact  was  accounted  for 
by  modifying  the  control-surface  input  coefficients  of  the  rigid-body 
equations  of  motion  to  include  static  aeroelastic  effects.  The  direc- 
tional flight-control  system  used  for  the  MIRA  simulation  includes  a 
rudder-flexibility  gain  to  account  for  this  phenomenon  (see  Section  3.8). 
However,  similar  information  is  not  available  for  6 , and  the  unmodified 

a 

data  are  used  as  presented  in  Section  3.4. 

The  final  equations  defining  the  three  lateral-directional  modes 
included  in  the  MIRA  simulation  are  given  in  Eq.  (3-100) , (3-101) , and 
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(3-102).  The  6^  and  6^  terms  present  in  these  equations  are  obtained 
in  the  same  fashion  as  the  5^  term  of  the  longitudinal  aeroelastic 
equations . 

' ^ = -0.831ii^  - 1566. 176n^ 


- ] -48.867p^gj^  - 50.649r„„T  + 0.0669n^  + 0.0755n- 
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- 7.6706  + 5.3266  - 1.420V, 

a r A 
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(3-102) 


3 . 8 Flight-Control  Systen. 

The  flight-control  system  implemented  in  the  MIRA  simulation  was 

derived  from  that  developed  for  the  Survivable  Flight  Control  System 
(2 ) 

(SFCS)  Program.  Although  the  goals  of  the  SF'CS  program  were  signi- 

ficantly different  from  those  of  the  MIRA  program,  the  flight-control 
system  developed  is  suitable  for  the  MIRA  investigation,  and  was  selected 
as  a representative  system.  The  flight-control  system  developed  for  the 
SFCS  program  is  an  adaptive  system  designed  for  use  over  the  whole  flight 
envelope  of  the  vehicle.  The  adaptive  portion  of  this  flight-control 
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system  has  not  been  programmed  in  the  MIRA  system,  but  reduced  to  a 
fixed-gain  system  suitable  for  operation  in  the  region  of  the  nominal 
flight  condition. 

Figures  3-5  and  3-6  depict  the  flight-control  system  programmed 
in  the  MIRA  simulation.  The  longitudinal  f light-control  system  was 
designed  to  provide  good  vehicle  handling  qualities,  and  used  a combina- 
tion of  An^  and  q feedback  to  achieve  this  goal.  The  lateral- 
directional  system  is  fairly  conventional  with  roll-rate  feedbacK  to 
the  ailerons,  and  washed-out  yaw-rate  and  lateral-acceleration  feedback 
to  the  rudder.  A roll-to-yaw  crossfeed  network  is  also  included  for 
turn  coordination.  In  addition,  it  is  possible  to  select  an  open-loop 
backup  flight-control  system,  should  it  be  desirable  to  investigate  this 
mode  of  operation.  The  flight-control  system  includes  structural-mode 
filters.  A second-order  notch  filter  centered  at  84  radians/second , 
followed  by  a first-order  lag  with  a corner  frequency  of  120  radians/ 
second,  is  included  in  the  longitudinal  channel.  Additional  phase  com- 
pensation is  also  needed.  The  roll  channel  includes  a second-order 
notch  filter  centered  at  40  radians/second.  In  the  directional  channel, 
the  only  structural-mode  filter  included  is  a first-order  lag  with  a 
corner  frequency  of  40  radians/second  on  the  lateral-accelerometer 
signal . 

The  longitudinal  SFCS  includes  a model  of  a special-purpose 
secondary  activator  which  has  been  deemed  to  be  overly  complex  for  the 
MIRA  task.  In  order  to  achieve  some  simplification  of  the  system,  it 
has  been  replaced  by  a second-order  model  which  is  identical  to  those 
included  in  the  roll  and  yaw  channels.  The  stick  force  gain  was  also 
modified  because  of  this  change  to  produce  a more  desirable  stick  force 
per  g.  In  all  three  channels,  the  power  cylinder  is  represented  by  a 
first-order  lag.  The  three  pilot-input  activating  signals  are  a center 
stick  for  longitudinal  control  with  a maximum  amplitude  of  50  pounds,  a 
center  stick  for  roll  control  with  limits  of  tl2  pounds,  and  the  rudder 
pedals  with  a maximum  force  input  of  150  pounds.  A feel  system  was  not 
mechanized,  thus  breakout,  nonlinearities,  etc.,  are  not  present  in  the 
MIRA  simulation. 

The  position  limits  of  the  aerodynamic  control  surfaces  are  also 
indicated  in  Figures  3-5  and  3-6.  They  are  +8  degrees  and  -20  degrees 
for  the  stabilator,  ±30  degrees  for  the  combined  aileron-spoiler  sur- 
face, and  ±30  degrees  for  the  rudder.  The  surfaces  are  also  rate- 
limited;  ±24  degrees/second  for  the  stabilator,  ±193  degrees/second  for 
the  aileron,  and  ±34  degrees/second  for  the  rudder. 
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Figure  3-5.  Longitudinal  flight-control  system. 


Figure  3-6.  Lateral-directional  flight-control  system. 


The  effect  of  rudder  flexibility  and  hinge-moment  limiting  is 
also  included  in  the  block  diagram  of  the  directional  system.  This 
phenomenon  is  modeled  by  a gain  less  than  1,  and  a position  limit  less 
than  or  equal  to  the  physical  limit  of  the  rudder  surface.  The  net 
effect  of  this  is  to  reduce  the  rudder  input  to  the  aircraft  aerodynamic 
and  structural-mode  equations.  At  the  nominal  flight  condition,  the 
rudder-flexibility  gain  is  0.86  and  the  hinge-moment  limit  is  ±15  degrees, 
rather  than  the  •:  30-degree  actual  position  limit. 

Separate  control  channels  are  mechanized  for  the  right  and  left 
engines.  A block  diagram  of  the  right  throttle  system  is  given  in 
Figure  3-7.  The  control  input  is  throttle  position,  which  is  limited 
to  values  between  0 and  100  degrees.  A gain  of  350  pounds/degree 
follows.  The  thrust  of  the  right  and  left  engines  is  limited  between 
0 and  35,000  pounds  in  each  case.  A first-order  lag  with  a time  con- 
stant of  0.5  second  follows  to  simulate  the  thrust  buildup  due  to  engine 
lags  following  a throttle-command  input. 


Figure  3-7.  Block  diagram  of  right-engine  control  s'^'tion. 

3 . 9 Simulation  Inputs 

The  ability  to  input  profile  time  histories  of  up  to  '.our  signals 
has  been  programmed  into  the  MIRA  simulation.  These  signals  can  be  any 
system  inputs — sucn  as  the  pilot  inputs,  the  gusts,  and  winds,  etc. 

This  capability  is  achieved  by  generating  a series  of  up  to  10  points 
in  time,  and  the  amplitudes  of  the  input  variables  at  the  corresponding 
points  in  time.  The  computer  program  linearly  interpolates  the  ampli- 
tude of  the  input  variables  between  the  time  points  to  achieve  the 
desired  profiles.  With  this  capability,  it  is  possible  to  apply  different 
kinds  and  combinations  of  inputs  to  the  simulation,  and  approximate  those 
present  in  actual  flight. 
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3.10 


Turbulence  Simulation 


A fairly  conventional  approach  to  the  simulation  of  turbulence  is 
used  in  the  MIRA  simulation.  The  technique  used  is  to  linearly  filter 
zero  mean  Gaussian  noise  to  generate  gust  signals  with  a Dryden  spec- 
trum. These  signals  are  easy  to  generate,  and  show  a reasonable 
agreement  with  measured  spectra.  Inherent  in  the  generation  of  gust 

signals  in  this  fashion  are  the  commonly  made  assumptions  that  the  gusts 
are  stationary  (the  statistical  properties  are  independent  of  time) , 
and  homogeneous  (the  statistical  properties  of  the  turbulence  field  are 
not  changed  by  a translation  of  the  coordinate  system  used  to  describe 
the  turbulence) . 


Since  a digital  program  is  available  which  generates  only  two 
independent  Gaussian  signals,  it  was  decided  to  generate  only 
and  Wg  in  the  MIRA  simulation.  The  u^  signal  has  been  dropped  since 
it  has  the  least  important  effect  of  the  three  on  the  aircraft  responses. 
The  side  and  vertical  gust  spectra  are  generated  by  passing  the  indepen- 
dent Gaussian  signals  through  filters  with  the  following  transfer  func- 
tions 


is  presently  fixed  at  the  trim  air  speed  for  the  MIRA  simulation. 


The  constant  is  needed  because  of  the  generation  of  the  gust  spectra 
by  digital  methods.  A gain  of  must  be  included  so  that  the 

This  assumes  that  the  variance 


For  the  nominal  flight  con- 
are  fixed  at  1,750  feet. 


variance  of  the  gust  signals  is  correct, 
of  the  Gaussian  generator  signal  is  unity, 
dition,  the  turbulence  scales,  and  L^, 

In  order  to  generate  stationary  random-gust  signals,  it  is  neces- 
sary to  independently  operate  the  gust-generation  scheme  before  applying 
the  signals  to  the  vehicle.  This  allows  the  system  transients  due  to  the 
filter  to  die  out,  and  is  accomplished  in  the  MIRA  simulation  by  switch- 
ing to  a standby  mode  prior  to  using  the  gust  signals  for  the  first 
time. 
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Several  statistical  analyses  subroutines  are  also  included  in  the 
MIRA  simulation.  It  is  possible  to  compute  the  mean,  standard  deviation 
and  mean-square  value  of  desired  signals.  Histograms  of  system  vari- 
ables can  also  be  generated. 


3.11  Sensor  Models 

In  Its  simplest  and  most  basic  terms,  the  MIRA  concept  involves 
the  use  of  aircraft  sensors  for  botli  navigation  and  flight-control 
purposes.  This  is  perhaps  the  most  important  and  unique  capability 
incorporated  into  the  MIRA  simulation.  The  sensor  signals  generated  in 
this  simulation  take  into  account  tlie  effect  of  location  on  the  accel- 
eration measurements,  gyro  inclination  effects,  and  aeroelasticity . 
Flight-control-system  sensor  models  arc  available  to  generate  signals 
for  flight-control  purposes  along  with  navigation-sensor  models  whose 
outputs  can  be  used  for  navigation,  as  well  as  flight-control,  purposes. 
At  present,  the  navigation  sensors  modeled  include  laser  gyros  and 
quartz  flexure  accelerometers.  This  section  contains  a discussion  of 
the  sensor  effects  and  models  incorporated  into  the  MIRA  simulation. 

The  rigid-body  accelerations  at  a particular  body-station  loca- 
tion relative  to  the  eg  are  given  by  the  following  matrix  equation. 

( [V]  + [w^J  [V]  - i[F  1 ) 

, ) 1 1 j grav  I 

[nl  = ^ I ; (3-105) 

•’0  1+  [;ox)[':l  + [i*ix  ] [wx  1 1 1 1 j 

Expansion  of  this  equation  to  allow  for  the  i2idependent  location  of  the 
three  accelerometers  leads  to 


+ g sin  0 - X (r^  + q^)  + y (pg)  + z (rp) 
0 L X XX. 
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The  flight-control  signal  uses  the  incremental  normal  acceleration  from 
trim,  and  the  initial  value  is  subtracted  from  n^  to  form  the  signal 
An^  for  this  purpose. 


The  capability  to  independently  locate  lateral  and  normal  accel- 
erometers for  flight  control  and  all  three  accelerometers  for  navigation 
has  been  programmed  into  the  MIRA  simulation.  In  practice,  it  is 
desired  to  locate  the  flight-control  normal  accelerometer  forward  of 
the  eg  so  that  the  stabilizing  effect  of  the  pitch-acceleration  compo- 
nent can  be  used  to  increase  the  short  period  damping.  The  location  of 
the  accelerometers  is  also  dictated  by  a desire  to  not  measure  unwanted 
accelerations.  For  example,  if  the  lateral  accelerometer  is  located  a 
significant  vertical  distance  from  the  waterline,  it  will  sense  unwanted 
roll  accelerations  in  addition  to  the  true  side  accelerations. 

The  rigid-body  rate-gyro  signals  are  not  affected  by  position, 
as  in  the  case  of  the  accelerometer  signals.  However,  the  roll-rate 
and  yaw-rate  gyros  sense  pitch,  roll,  and  yaw  rates  due  to  gyro  inclina- 
tion effects.  To  minimize  cross-axes  sensing,  the  roll-  and  yaw-rate 
gyro  sensitive  axes  are  located  parallel  to  and  perpendicular  to  the 
average  longitudinal  principal  axis,  respectively.  The  longitudinal 
principal  axis  is  oriented  an  average  of  1.5  degrees  below  the  aircraft 
waterline  about  which  the  body-axes  equations  are  written.  In  order 

to  obtain  the  roll-rate  and  yaw-rate  signals  about  the  principal  axis, 
it  is  necessary  to  transform  p and  r through  a rotation  of  -1.5  degrees, 
thus 

Pg  = p cos  - r sin  (3-109) 

r,,  = p sin  a„  + r cos  a_  (3-110) 

L>  o o 

Bending  effects  are  incorporated  into  the  sensor  signals  by  the 
addition  of  appropriate  terms  to  the  rigid-body  signals.  Acceleration 
signals  are  corrupted  by  the  product  of  the  displacement  due  to  the 
mode  and  the  second  derivative  of  the  generalized  mode  shape.  Acceler- 
ometers are  optimally  located  at  the  nodes  of  the  modes  to  eliminate 
the  effects  of  aeroelasticity . On  the  other  hand,  rate  gyros  are 
affected  by  the  product  of  the  slope  of  the  mode  shape  and  the  rate 
of  the  generalized  mode  shape.  Their  optimum  location  is  at  the  anti- 
nodes of  the  mode  shapes.  In  practice,  it  is  generally  impossible  to 
locate  sensors  at  their  optimum  locations  because  of  physical  constraints 
and  the  need  to  satisfy  the  constraints  of  each  mode  simultaneously. 

Based  on  the  data  presented  in  References  2 and  3,  the  following 
equations  result  (for  the  sensor  signals  which  include  the  effects  of 
bending) . 


44 


r ^ 


(3-111)  ; 

(3-112)  I 

(3-113) 

i 

(3-114) 

(3-115) 

The  coefficients  of  the  bending  terms  in  the  previous  equations  are  a 
function  of  the  location  of  the  sensors.  The  yaw-rate-gyro  signal  and 
the  lateral-accelerometer  signal  are  affected  only  by  n,,  the  fuselage 

D 

first  lateral-bending  mode,  and  not  by  the  fuselage  first  torsion  mode 
and  the  wing  first  asymmetric-bending  mode.  The  roll-rate  signal  is 
affected  by  all  three  lateral-directional  modes.  Each  of  the  three 
longitudinal  modes  also  corrupts  the  pitch-rate-gyro  and  normal- 
acceleration  signals.  The  longitudinal-accelerometer  signal  is  not 
affected  by  any  mode  signals  for  reasons  discussed  previously. 

The  flight-control-system  sensor  models  incorporated  into  the 
MIRA  simulation  are  those  used  for  the  SFCS  simulation  and  analysis. 

In  the  case  of  the  fight-control-system  accelerometers,  a first-order 
lag  represents  the  dynamics  of  the  sensors.  The  corner  frequency  of 
the  lag  is  200  radians/second  for  the  normal  accelerometer,  and 
220  radians/second  for  the  lateral  accelerometer.  Second-order  systems 
are  used  to  represent  the  flight-control-system  rate  gyros.  The 
frequency  of  these  systems  is  150  radians/second,  and  the  damping  ratio 
is  0.6.  These  sensor  models  operate  on  the  signals  given  by  Eq.  (3-111) 
through  (3-115)  to  produce  the  input  signals  to  the  flight-control 
system. 

Limits  are  imposed  on  the  rates  which  the  flight-control-system 
gyros  sense.  The  limits  are  ±200  degrees/second  for  the  roll-rate  gyro, 

±50  degrees/second  for  the  pitch-rate  gyro,  and  ±150  degrees/second  for 
the  yaw-rate  gyro,^^^  The  flight-control  lateral  accelerometer  is 
limited  to  ±2  g,  while  the  normal  accelerometer  is  limited  to  10  g up 
and  4 g down. 
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More  sophisticated  models  are  incorporated  into  the  MIRA  simula- 
tion for  the  navigation  gyros  and  accelerometers.  Error  models  for 
three  orthogonal  ring-laser  gyros  are  programmed.  Each  has  a maximum 
rate  of  400  degrees/second  and  a maximum  pulse  rate  of  916,000  pulses/ 
second. 


Figure  3-8  is  a block  diagram  of  the  roll-rate  gyro,  which  is 
typical  of  the  other  two  as  well.  The  first  operation  undertaken  in 
the  model  is  the  integration  of  the  roll-rate-gyro  signal,  including 
bending  effects.  This  is  done  by  rectangular  integration,  i.e. 
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Figure  3-8.  Block  diagram  of  ring-laser  roll-rate  gyro. 


The  resultant  signal  represents  the  output  of  a perfect  gyro.  Next, 
misalignment  errors  are  introduced  with  the  addition  of  terms  propor- 
tional to  the  integrated  pitch-rate  and  yaw-rate  signals. 
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(3-117) 


PgIM  ^GI  ^l^GI  '^2’^GI 


Multiplication  by  a scale  factor  converts  the  signal  to  pulses  and  a 
bias-error  signal  is  added  in.  The  computational  algorithm  is 


■ GCT 


= ^^p^<Pgim> 


+ Ab 


(3-118) 


It  has  been  determined  that  it  is  not  necessary  to  model  the  random-wal)c 
effects  which  characterize  laser-gyro  sensors  for  this  phase  of  the  MIRA 
project  because  it  is  a long-term  phenomenon.  When  it  is  decided  to 
inco-porate  this  effect,  it  will  be  summed  in  at  this  juncture  in  the 
model.  Quantization  effects  are  represented  by  adding  the  fractional 
part  of  the  pulse  count  remaining  from  the  previous  iteration  to  the 
present  pulse  count,  and  integerizing  the  sum.  The  result  is  called 
Pgcti"  navigation-system  signal  is  obtained  by  introducing  the 

effects  of  the  scale  factor  and  its  error,  i.e. 

^p JPgCTI  (3-119) 

z Pg  ' G/ 


An  option  is  included  in  the  simulation  to  obtain  the  error  between 
the  perfect  navigation  sensor  and  that  corrupted  by  all  of  the  errors 
discussed.  Tlie  flight-control  roll-rate  signal  is  obtained  from  this 

sensor  model  by  dividing  AOg  by  the  computation  time  increment. 

z 

Three  orthogonal  quartz-flexure-type  accelerometers,  are  modeled 
in  the  MIRA  simulation.  Each  instrument  can  sense  a maximum  linear 
acceleration  of  10  g.  A maximum  pulse  rate  of  320,000  pulses/second, 
and  correspondingly  relatively  fine  quantization  is  achieved  through 
the  use  of  the  Sigma-Delta  Modulation  technique.  Figure  3-9  is  a 
bloc)t  diagram  of  the  lateral-accelerometer  model,  which  is  typical  of 
the  other  two  instruments.  First,  the  lateral  acceleration  including 
bending  is  rectangularly  integrated  to  produce  a perfect  accelerom- 
eter signal,  i.e. 
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Misalignment  errors  are  modeled  by  adding  terms  proportional  to  the 
integrated  normal-  and  longitudinal-accelerometer  signals  to  the  inte- 
grated lateral-accelerometer  signal.  That  is 
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Figure  3-9.  Block  diagram  of  quartz-flexure  lateral  accelerometer. 
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The  resultant  signal  is  then  passed  through  a first-order  lag  to  repre- 
sent delays  in  the  instrument 
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Crosscoupling  effects  are  introduced  by  adding  a constant  times  the  pro- 
duct of  the  integrated  lagged  longitudinal-accelerometer  signal,  including 
misalignment  effects,  and  the  latest  lateral-accelerometer  signal  to  the 
lateral-accelerometer  signal.  Nonlinearities  are  taken  into  account  by 
summing  in  a constant  times  the  square  of  the  integrated  lagged  lateral- 
acceleration  signal.  A term  representing  the  bias  error  is  also  intro- 
duced, and  the  resultant  signal  converted  to  a pulse  count  by  multiply- 
ing by  a scale  factor.  The  algorithm  representing  these  operations  is 
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Quantization  effects  are  handled  by  summing  in  the  fractional  part  of 
the  pulse  count  remaining  from  the  previous  iteration  to  the  present 
pulse  count  and  integerizing  the  sum,  to  produce  a signal  designated  as 
^CTI‘  navigation  signal  is  obtained  by  reintroducing  the  effects 

of  the  scale  factor  along  with  its  error,  i.e. 


= SP^(l  + ^n  ji'cTI  (3-124) 
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An  option  is  included  in  the  simulation  to  obtain  the  error  between  the 
perfect  navigation  signal  and  the  one  corrupted  by  all  of  the  errors 
discussed.  The  acceleration  signal  for  flight  control  is  obtained  by 
dividing  the  navigation  acceleration  signal  by  At. 


3.12  Navigation  System  Model 

Presented  in  this  section  is  the  algorithm  programmed  in  the  MIRA 
simulation  representing  a strapdown  Local  Vertical  Wander  Azimuth  (LVWA) 
navigation  system.  Reference  5 should  be  consulted  for  a more  detailed 
description  of  the  derivation  of  the  algorithm  and  the  assumptions 
involved  in  its  generation.  The  navigation  model  can  be  divided  into 
two  major  submodules;  a high-speed  loop  computed  at  50  hertz,  and  a 
low-speed  loop  computed  at  5 hertz. 
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The  high-speed  navigation  loop  involves  the  generation  of  the 
body-to-computational  or  LVWA  transformation  matrix,  the  transformation 
of  the  navigation  sensor  velocities  from  the  body  frame  to  the  computa- 
tional frame  and  their  subsequent  accumulation  over  the  low-speed  com- 
putation time  interval.  A first-,  second-  or  third-order  quaternion 
update  can  be  used  to  generate  the  body-to-computational-f rame  trans- 
formation matrix.  As  recommended  in  Reference  5,  the  generation  of  the 
body-to-computational-f rame  transformation  matrix  and  the  transformation 
of  the  velocity  signals  are  staggered  in  time.  The  body-to-computational- 
frame  transformation  matrix  is  computed  midway  between  the  times  at  which 
the  velocity  transformation  occurs.  The  attitude  angles  are  extracted 
from  the  body-to-computational-f rame  transformation  matrix  for  compari- 
son with  the  corresponding  true  angles  generated  in  the  vehicle  portion 
of  the  simulation. 

The  function  of  the  low-speed  portion  of  the  navigation-system 
model  is  the  integration  of  the  fundamental  equation  of  inertial  navi- 
gation for  a navigator  employing  an  earth-relative  computational  frame. 
This  equation  is 

[Vj^c'  = (3-125) 

All  signals  are  defined  in  the  computational  frame.  A second  integra- 
tion follows  to  obtain  the  earth-relative  position  of  the  vehicle. 

Since  the  earth- relative  position  is  expressed  in  angular  coordinates 
(latitude,  longitude,  and  wander  angle) , it  is  necessary  to  convert  the 
level  components  of  linear  velocity  from  the  first  integration  into 
angular  velocity  by  dividing  by  the  radius  of  curvature.  The  result  is 
the  angular  velocities  of  the  computational  frame  with  respect  to  the 
earth.  The  computation-f rame-to-earth  transformation  matrix  is  then 
obtained  by  integrating  the  following  matrix  differential  equation 

[TceI  - (3-126) 

and  latitude,  longitude,  and  wander  angle  extracted  from  via 

inverse  trigonometric  functions. 

Altitude  is  obtained  by  integrating  the  vertical  component  of  the 
earth-relative  velocity.  Third-order  damping  is  included  to  stabilize 
the  system.  The  navigation  algorithm  of  Reference  5 has  been  modified 
to  reflect  the  fact  that  a spherical  earth  is  used  in  the  vehicle 
equations  of  motion.  This  step  involved  the  use  of  a gravity  model 
different  from  that  given  in  the  reference,  and  the  use  of  a constant 
earth  radius  to  compute  the  angular  velocities  about  level  axes. 
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A transformation  of  the  navigation-sensor  signals  is  required 
before  use  in  the  navigation  algorithm.  This  is  due  to  the  use  of  a 
different  coordinate  frame  for  the  definition  of  the  navigation  algo- 
rithm from  that  used  to  define  the  vehicle  equations  of  motion.  The 
X,  y,  z coordinates  of  the  navigation  system  are  defined  in  an  up,  east, 
north  sense — -while  those  of  the  vehicle  are  in  a nominal  north,  east, 
down  sense.  Thus,  the  required  transformation  from  the  body  frame  to 
the  navigation  frame  is 
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'B 


= Yr 


(3-127) 

(3-128) 

(3-129) 


The  airframe  equations  with  bending  and  the  sensor  equations  are 
solved  at  400  hertz,  while  the  high-speed  navigation  loop  is  solved  at 
50  hertz.  It  is  necessary  to  sum  the  sensor  outputs  over  the  navigation 
system  computation  interval  to  account  for  this  difference  in  the  fre- 
quency of  computation. 

The  equations  of  the  high-speed  navigation  loop  are  now  considered 
in  detail.  The  integral  of  the  angular  rate  of  the  body  frame  with 
respect  to  the  computational  frame  is  approximated  by  tailing  the  inte- 
grated angular  rates  of  the  body  frame  with  respect  to  the  inertial 
reference  (the  outputs  of  the  navigation  sensors) , and  subtracting  the 
integrated  angular  rates  of  the  computational  frame  with  respect  to  the 
inertial  reference  in  body-frame  coordinates.  The  integrated  angular 
rates  of  the  computational  frame  with  respect  to  the  inertial  reference 
are  obtained  from  the  low-speed  navigation  loop,  and  are  the  sum  of  the 
inertial  earth-rate  signals  and  the  computation-frame  signals  with 
respect  to  earth  times  the  high-speed  computation  time.  These  signals 
must  then  be  transformed  from  the  computational  frame  to  the  body  frame. 
This  formulation  assumes  that  the  angular  accelerations  of  the  computa- 
tion frame  with  respect  to  the  inertial  reference  are  zero.  In  the 
simulation,  the  following  equations  are  solved. 
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(3-130) 
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Next,  the  parameters  of  the  quaternion  update  are  computed  using 
the  integrated  angular  rates  of  the  body  frame  with  respect  to  the  com- 
putational frame  just  obtained.  Either  a first-,  second-,  or  third- 
order  update  can  be  used,  and  the  appropriate  set  of  the  following  param- 
eters computed. 


First-order  update 


Pq  = 1 (3-132) 
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Second-order  update 
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Third-order  update 

Po  = 

1 2 
1 - |(A9^) 

(3-141) 
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(3-142) 
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(3-144) 

The  subscript  n in  the  previous  ecjuations  denotes  that  the  values  from 
the  present  iteration  are  used,  while  the  variables  witn  a subscript 
r-1  are  those  from  the  previous  iteration.  Wiion  subscripts  are  not 
used,  the  present  values  are  used.  The  updated  quarternion  parameters 
are  then  given  by 
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Once  every  second,  the  q terms  are  normalized  by  multiplying  each  one 
by 
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The  initial  values  of  the  q terms  are  obtained  from  the  expressions 
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where  4>j./  9^. « and  are  the  trim  values  of  the  roll,  pitch,  and  yaw 

angles;  and  the  initial  wander  angle  is  assumed  to  be  zero. 


With  respect  to  the  q terms,  the  transformation  from  body  to  com- 
putational coordinates  is  given  by 
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whereas,  in  terms  of  the  Euler  angles  this  matrix  is 
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'^AN‘  navigation  Euler  angles  can  then  be  obtained  by 


the  following  formulae,  for  comparison  with  the  vehicle's  Euler  angles. 
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A transformation  of  the  integrated  linear  accelerations  from  the 
body  frame  to  the  computational  frame  is  performed  next,  using  tne  body- 
to-computaticnal-f rame  transformation  matrix  just  defined.  The  appro- 
priate equations  are 
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Then,  these  signals  are  summed  over  the  low-speed  navigation-system  com- 
putation interval  before  being  fed  to  the  low-speed  portion  of  the 
navigation  system,  i.o. 
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This  completes  the  discussion  of  the  high-speed  portion  of  the 
navigation-system  model. 

A detailed  discussion  of  the  algorithms  of  the  low-speed  portion 
of  the  navigation  model  follov;s.  The  first  step  is  the  updating  of  the 


computational-frame  velocities.  The 
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The  lii^r  and  are 

the 

elements  of  the 

matrix  ] 

[V„_]  • At._,  which  results  from  the  fundamental  equation  of  the  navi- 

NC,  LiU 

gation  system.  The  initial  values  of  , and  are  obtained 

by  transforming  the  initial  values  of  u , , v , , and  w , to  the  com- 
■*  rei  rel  rel 

putational  frame. 

Damping  of  the  vertical  velocity  is  accomplished  by  recalculating 
the  vertical-vector  component  according  to  the  equation 
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Where  6h  = h - h„  is  the  altitude  error,  and  6a  is  defined  later  in  the 
N 

discussion  of  the  low-speed  navigation  system  when  the  altitude  is 
calculated . 

The  calculation  of  the  angular  velocities  about  the  level  axes 
occurs  next.  For  a spherical  earth,  the  appropriate  equations  are 
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An  updating  of  the  computational-f rame-to-earth  direction  cosine 
matrix  follows.  This  is  accomplished  by  solving  the  equation 
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The  matrix  [T__]  is  given,  in  terms  of  lattitude,  longitude,  and  wander 
angle  by 
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and  is  initialized  at  w^  = 0.0,  X^,  and  o^.  Every  30  seconds,  the 
direction  cosine  matrix  is  orthonormalized  according  to  the  formula 

i, 


r T 

L ortho 


/^Ce]  “ 2['^CeJ  )[^CE_  7ceJ  ” l/o.  I (3-170) 


The  computation  of  the  earth-rate  signals  occurs  next  in  the 
navigation  model,  using  the  equations 
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The  gyro-torquing  signals  required  for  the  generation  of  the  integrated 
angular  rates  of  the  body  frame,  with  respect  to  the  computational  frame 
in  the  high-speed-navigation  loop,  can  now  be  obtained.  These  signals 
are  generated  from  the  equation 

[wci]  = (3-172) 

or  more  explicitly 
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The  integrated  angular-rate  terms,  crosscoupled  with  the 
computational-frame  velocity  terms  present  in  the  velocity  update 
equation  (u^,  and  , are  now  calculated  using  the  expression 
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expressions  for  the  elements  of  this  matrix  are 
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Computation  of  the  navigation-system  altitude  is  accomplished  by 
integrating  the  vertical-velocity  component.  Including  the  effects  of 
vertical  damping,  the  appropriate  equations  are 
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The  altitude  error  is  also  updated  in  this  portion  of  the  model  for  use 
in  the  next  cycle  of  the  damping  of  the  vertical-velocity  term. 
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Next,  the  gravity  terms  included  in  the  velocity-update  equation 
are  calculated  from  expressions  which  include  the  central  and  centri- 
fugal forces  because  of  the  spherical-earth  approximation.  The  appro- 
priate equations  are 
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Latitude,  longitude,  and  wander  angle  can  be  extracted  from  the 
elements  of  the  direction  cosine  matrix,  viz 
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Last  of  all,  the  sine  and  cosine  of  the  wander  angle  are  used  to 
transform  the  navigation  velocities  in  the  computational  frame  to  the 
local-vertical  north-pointing  frame,  according  to  the  formulae 


(d) 


(3-188) 
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NC  AN 
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To  aid  in  the  evaluation  of  different  MIRA  concepts,  error  signals 
between  vehicle  and  navigation-system  variables  are  generated  in  the 
simulation.  The  errors  generated  are  those  of  the  Euler  angles,  lati- 
tude, longitude,  wander  angle,  altitude,  and  the  earth-relative 
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velocities.  This  completes  the  discussion  of  the  computations  performed 
in  the  LVWA  navigation-system  mode  of  the  MIRA  simulation. 

3.13  A General  Description  of  the  Navigation  and  Attitude-Error 
Propagation  Model 

This  section  contains  a general  description  of  an  algorithm  devel- 
oped to  generate  the  error  in: 


(1) 

Position 

— latitude,  longitude,  altitude 

(2) 

Velocity 

— north , east , 

vertical . 

(3) 

Attitude 

— roll,  pitch. 

heading . 

for  a strapdown  system.  It  is  assumed  that  the  incremental  velocity 
and  incremental  angular  errors  are  available  from  cascaded  inertial- 
sensor models  and  compensation,  or  inertial-sensor  error  models, 
together  with  the  usual  flight-profile  data.  This  teclinique  (described 
in  detail  in  Reference  8)  will  replace  that  presently  incorporated  into 
the  MIRA  simulation,  which  includes  the  strapdown-navigation  model  and 
subsequent  differencing  of  its  outputs  with  the  aircraft  equations-of- 
motion  model  in  order  to  generate  the  errors  discussed  in  Section  3.12. 
The  approach  currently  implemented  in  the  MIRA  simulation  more  accu- 
rately reflects  the  expected  performance  of  an  actual  INS,  since  it 
includes  the  effects  of  computational  errors  in  the  overall  navigator 
performance.  However,  it  also  introduces  more  stringent  computational 
requirements  than  are  necessary  to  evaluate  the  effects  of  the  flight 
environment  on  the  inertial-sensor  performance  and  the  overall  system 
performance.  In  the  proposed  program,  the  only  differencing  remaining 
is  that  between  the  simulated  inertial  sensors  and  the  "ideal"  sensors. 
This  feature  renders  the  proposed  scheme  different  from  any  existing 
error-propagation  program. 

The  following  assumptions  are  made  in  the  development  of  the 
navigation-error  equations: 

(1)  The  earth  is  spherical. 

(2)  The  aircraft  equations-of-motion  model  provides  a nominal 
trajectory,  consisting  of  the  body-to-local-vertical  trans- 
formation, the  body  specific  forces,  the  body  angular  rates, 
the  latitude,  the  longitude,  and  the  altitude  of  the  vehicle, 
and  the  local  vertical  with  respect  to  earth  velocities. 

(3)  The  inertial-sensor  models  provide  the  errors  in  the  angular 
rates  and  specific  forces. 
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since  the  earth  model  used  in  both  the  equations  of  motion  and 
the  strapdown-navigation  model  is  spherical,  a space-stabilized  compu- 
tational frame  is  employed  for  the  navigation-error  equations.  This 
has  the  great  advantage  that  the  navigation  equations,  including  verti- 
cal damping,  are  perfectly  symmetrical  with  respect  to  each  of  the  three 
axes  of  the  position  and  velocity  variables.  Therefore,  only  one  axis 
of  the  three  need  be  considered  in  developing  the  error  equations,  and 
the  results  extended  to  the  remaining  two  axes  by  symmetry.  The  result 
is  that  the  state  errors  are  arranged  in  a 2x3-matrix  form  rather  than 
as  a 6xl-column  vector,  and  the  state  transition  matrix  becomes  2^2, 
instead  of  6x6.  Similarly,  the  attitude-error  equations  become  a single 
{3x3)x(3xl)  matrix-vector  equation. 

The  attitude  computations  for  a strapdown  system  consist  merely 
of  extracting  the  three  Euler  angles  from  the  body-to-local-vertical- 
north  (LVN)  transformation  via  appropriate  inverse  trigonometric 
functions . 

Attitude  errors  arise  from  two  different,  but  not  independent, 
sources.  First,  any  errors  in  initial  alignment  or  gyro-drift  compen- 
sation will  produce  angular  errors  in  the  knowledge  of  the  motion  of 
the  body  frame.  Second,  accelerometer-  and  gyro-compensation,  initial- 
alignment,  and  mechanization  errors  all  integrate  into  latitude  and 
longitude  errors.  These  are  angular  errors  in  the  knowledge  of  the 
LVN  frame.  Since  the  computed  attitude  matrix  is  a function  of  both 
where  the  computer  "thinks"  the  body  frame  is  and  where  it  "thinks" 
the  LVN  frame  is,  the  attitude-error  vector  is  a function  of  these  two 
factors,  as  well  as  the  attitude  angles.  Naturally,  the  attitude-error 
computations  are  singular  for  a pitch  angle  of  *90  degrees.  In  deriving 
the  attitude-error  equations,  the  position  errors  must  first  be  inter- 
preted as  errors  in  the  transformation  from  the  inertial  frame  to  the 
computed  LVN  frame.  Then,  the  errors  in  the  LVN  transformation  must 
be  interpreted  as  attitude  errors.  Finally,  the  matrix  of  the  result- 
ing system  of  equations  must  be  inverted  to  obtain  the  attitude  errors 
explicitly . 

The  main  savings  resulting  from  the  error-propagation  approach 
are  the  elimination  of  the  strapdown  "velocity/attitude"  algorithm,  and 
the  elimination  of  the  attitude  extraction.  In  addition,  the  navigation- 
like computations  for  a strapdown  system  are  more  efficiently  performed 
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in  the  inertial  frame.  Finally,  computational-accuracy  requirements 
are  reduced,  since  errors  are  the  only  forcing  functions. 

In  addition  to  the  efficiency  resulting  from  this  approach,  it 
isolates  the  effects  of  sensor  errors  from  all  other  navigation  effects 
(truncation,  roundoff,  model  limitations),  and  thus  is  ideally  suited 
for  the  MIRA  application. 
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SECTION  4 


DEMONSTRATION  OF  MIRA  SIMULATION  CAPABILITY  AND 
COMPARISON  OF  RESULTS 


4 . 1 Introduction 

This  section  deals  with  a demonstration  of  the  simulation  capa- 
bility developed  to  evaluate  the  MIRA  concept.  First,  a description 
of  the  sensor  configurations  selected  for  evaluation  is  given.  Then, 
associated  data  (such  as  the  structural-mode  sensor  coefficients,  the 
pilot  inputs  to  the  simulation,  the  turbulence  levels,  the  navigation- 
sensor parameters,  etc.)  are  discussed.  A general  comparison  of  the 
results  obtained  for  the  various  cases  follows. 

The  data  contained  in  this  section  should  not  be  viewed  as  repre- 
sentative of  a recommended  MIRA  configuration,  but  as  a demonstration 
of  the  capability  of  the  simulation.  Several  configurations  were  chosen 
which,  at  this  point,  appear  to  be  valid  candidates  for  MIRA  systems 
(i.e.,  both  dispersed  and  colocated  inertial  instruments,  and  both 
separate  and  common  navigation  and  flight-control  instruments)  to 
demonstrate  the  feasibility  of  using  the  simulation  for  a wide  class 
of  MIRA  systems.  In  general,  the  test  cases  included  in  this  report 
are  not  of  sufficient  duration  to  permit  the  drawing  of  substantial 
conclusions  concerning  the  navigation  performance  of  the  various  con- 
figurations. The  test  run  times  were  limited  to  allow  a wide  range  of 
cases  to  be  investigated  in  this  pha^e  with  the  non-real-time  simula- 
tion. Within  this  context,  however,  some  limited  conclusions  are  drawn 
about  the  MIRA  configurations  considered,  and  about  the  MIRA  approach 
in  general  as  it  impacts  navigation  and  flight-control  performance. 

The  importance  of  the  bending-mode  data  in  this  performance  evaluation, 
particularly  in  the  area  of  flight  control,  is  clearly  indicated  by  the 
data  to  be  presented  herein. 

4 . 2 Data  for  MIRA  Simulation  Runs 

Three  basic  sensor  configurations  have  been  selected  for  evalua- 
tion in  this  phase  of  the  project.  The  configuration  used  to  check  out 
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the  simulation,  as  well  as  to  demonstrate  the  MIRA  simulation  capability, 
had  sensors  dispersed  along  the  fuselage  at  locations  selected  for  the 
flight  test  of  the  Survivable  Flight  Control  System. This  is 
denoted  as  the  SFCS  configuration.  The  sensor  locations  for  this  con- 
figuration were  chosen  based  upon  space  available  on  the  test  aircraft 
to  mount  the  sensors,  and  a desire  to  minimize  the  effect  of  the  aero- 
elastic  modes  on  the  sensor  responses  given  these  physical  constraints. 
The  sensor  locations  for  this  configuration  are  given  in  Table  4-1.  As 
a point  of  reference,  the  eg  for  the  nominal-vehicle  configuration  is 
located  at  30.6  percent  c,  or  320  inches  behind  the  nose  of  the  aircraft 
along  the  waterline  (FS320),  and  27.65  inches  up  from  the  waterline  of 
the  vehicle  (WL27.65). 


Table  4-1.  SFCS  sensor  locations. 


Sensor 

Fuselage 
Station 
(FS) 
(in. ) 

VJater- 
Line 
(WL) 
(in. ) 

Butt- 
line 
(in. ) 

Direction  of 
Sensing  Axis 

Pitch-Rate  Gyro 

383 

6 

47.5  L 

Parallel  to  the  y axis, 
positive  aircraft  nose 
up 

Roll-Rate  Gyro 

77 

11.5 

00 

Forward  in  the  x-z  plane 
and  1.5  degrees  below 
waterline 

Yaw-Rate  Gyro 

383 

6 

47.5  R 

Down,  parallel  to  x-z 
plane  and  perpendicular 
to  the  roll-rate-gyro 
sensing  axis 

Lateral 

Accelerometer 

186.3 

19.8 

00 

Perpendicular  to  x-z 
plane,  positive  to 
pilot's  right 

Normal 

Accelerometer 

77 

22.5 

00 

Perpendicular  to  the 
waterline,  positive  down 

The  other  two  configurations  selected  for  evaluation  assume  that 
all  of  the  sensors  are  colocated.  One  sensor  package  was  chosen  to 
be  located  slightly  forward  of  the  eg  at  FS313  on  the  x body  axis,  and 
27.65  inches  up  from  the  waterline.  This  location  was  selected  because 
it  was  one  of  the  candidate  positions  for  locating  sensors  on  the  SFCS 
program,  and  presumably  there  is  some  available  space  there.  The  third 
sensor  location  selected  for  evaluation  in  this  phase  of  the  program 
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assumes  that  all  of  the  sensors  are  located  just  forward  of  the  pilot's 
station  at  FS77.  The  pilot's  station  is  located  at  about  FS119.  Relative 
to  the  waterline,  this  package  is  assumed  to  be  located  at  WL22.5.  For 
all  configurations,  the  roll  rate  gyro  sensing  axis  is  inclined  1.5°  below 
the  waterline.  The  yaw  rate  gyro  sensing  axis  is  perpendicular  to  the 
roll  rate  gyro  sensing  axis. 

This  selection  of  sensor  configurations  allows  many  different  cases 
to  be  evaluated.  It  is  possible  to  investigate  system  performance  using 
any  one  of  the  three  configurations  for  the  flight-control  sensors,  while 
simultaneously  using  any  one  of  the  configurations  for  the  navigation  sen- 
sors. The  use  of  common  sensors  for  both  navigation  and  flight-control  pur- 
poses can  also  be  investigated  for  each  of  the  three  selected  configurations. 
Furthermore,  perfect  sensors  can  be  selected,  and  the  level  of  turbulence 
varied. 

Table  4-2  contains  a list  of  the  configurations  evaluated  using 
the  MIRA  simulation.  The  quantity  of  data  gathered  from  these  runs  was 
quite  vast.  In  order  to  reasonably  limit  the  amount  of  data  contained 
herein,  while  at  the  same  time  demonstrating  the  capability  of  the  MIRA 
simulation,  results  are  presented  to  allow  a comparison  of  the  selected 
sensor  configurations  from  two  points  of  view.  In  the  first  instance, 
a basic-sensor  configuration  is  selected;  e.g.,  the  sensors  are  at  the 
SFCS  locations,  A comparison  is  then  made  for  three  different  condi- 
tions. The  first  is  an  ideal  situation  where  perfect  sensors  are  used 
for  navigation  and  flight-control,  and  the  vehicle  is  flying  in  a non- 
turbulent  environment.  In  the  second,  the  perfect-sensor  models  are 
replaced  by  realistic  models,  and  the  vehicle  is  flying  in  a turbulent 
environment.  For  the  third  condition,  tlie  navigation  sensors  are  used 
for  both  navigation  and  flight  control,  and  the  vehicle  is  flying  in 
gusts.  More  specifically,  the  basic  responses  of  the  vehicle  and  naviga- 
tion system  for  configurations  11,  12,  and  13  (see  Table  4-2)  are  com- 
pared, In  addition,  configurations  8,  9,  and  10  (see  Table  4-2)  are 
compared  with  regard  to  their  structural-mode  responses. 

A comparison  is  also  made  of  the  vehicle  responses  obtained  with 
a conventional  sensor  configuration  (in  which  the  flight-control  sensors 
are  located  at  the  SFCS  stations  on  the  vehicle,  and  the  navigation 
sensors  are  located  close  to  the  eg  of  the  vehicle  at  FS313).  These 
responses  are  compared  with  those  obtained  for  three  sensor  configura- 
tions using  common  sensors  for  both  navigation  and  flight-control  pur- 
poses. In  one  case,  the  sensors  are  at  FS313;  in  another,  forward  of 
the  pilot's  station  at  FS77;  and,  in  the  third,  at  the  SFCS  location. 
Configurations  3,  7,  10,  and  13  (see  Table  4-2)  are  the  specific  ones 
compared  in  this  report. 
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rms  value  of  10  ft/s. 


A fairly  simple  maneuver  is  used  for  the  purpose  of  demonstrating 
the  capability  of  the  MIRA  simulation.  The  specific  pilot  inputs  to 
the  simulation  are  shown  in  Figure  4-1.  A right  lateral-stick  pulse 
input  of  1.33  pounds  is  applied  to  the  vehicle  in  order  to  achieve  a 
banked  right  turn.  A bank  angle  of  approximately  32.5  degrees  results. 
Simultaneously,  a pitch-stick  step  input  of  0.915  pound  is  applied  to 
keep  the  nose  of  the  aircraft  up.  These  are  the  inputs  to  the  vehicle 
until  a heading  change  of  approximately  37.5  degrees  is  achieved,  at 
which  time,  the  pitch-stick  input  is  removed  and  a left  lateral-stick 
pulse  of  1.33  pounds  is  applied  to  achieve  wings-level  flight.  After 
a short  period  to  allow  the  system  transients  to  die  out,  a pitch-stick 
step  input  of  5 pounds  is  applied  to  pull  the  nose  of  the  aircraft  up. 
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The  coefficients  defining  the  effect  of  the  aeroelastic  inodes  of 
the  F-4  on  the  sensor  responses  were  determined  for  each  of  the  sensor 
configurations  discussed  previously.  Tables  4-3  and  4-4  contain  the 
pertinent  parameters,  along  with  the  distances  from  the  accelerometer 
locations  to  the  vehicle  eg.  These  data  were  obtained  from  References  2 
and  3,  but  no  information  is  given  about  the  location  of  the  longitudinal 
accelerometer,  since  it  is  not  needed  for  flight-control  purposes.  How- 
ever, for  the  SFCS  sensor  configuration,  it  is  assumed  that  the  longi- 
tudinal accelerometer  is  located  at  FS313,  so  that  each  of  the  three 
accelerometers  has  a different  physical  location. 

Table  4-3.  Aeroelastic  coefficients  of  acceleration  responses. 


^ mated  coefficient. 


Table  4-4.  Aeroelastic  coefficients  of  gyro  responses. 


Pitch-Rate  Gyro 

Sensor 

Configuration 

30/3nj^  (rad) 

36/3n2  (rad) 

36/3ri3  (rad) 

FS77 

SFCS (FS383) 
FS313 

-3.608(10”'') 

0.0 

-2.33  (10”'') 

2.239  (lO”^) 
-7.11(10"'*) 
-2.64  (10”^) 

_4  * 

-1.14(10  ) 

1.14  dO”^) 
0.0 

Roll-Rate  Gyro 

Sensor 

Configuration 

3<I>/3n^  (rad) 

3((i/3n5  (rad) 

3ii/3ri^  (rad) 

FS77 

SFCS (FS77) 
FS313 

-5.5  do"'') 
-5.5(10"^) 
-3.33(10”^) 

6.3  dO”'') 
6.3  (10"“') 
5.333  (10”“') 

2.9(10”^) 
2.9(10”^) 
2.267  dO”^) 

Yaw-Rate  Gyro 

Sensor 

Configuration 

3i);/3ng  (rad) 

FS77 

-1.4(10  ■") 

SFCS (FS383) 

6.0(10”^) 

FS313 

-3.5  dO”^) 

* 

Estimated  coefficient. 


The  normal  accelerometer  is  located  at  either  FS77  or  FS313,  and 
the  coefficient  of  the  first  longitudinal-structural  mode  changes  sign 
between  these  two  locations,  with  its  magnitude  at  FS77  about  three 
times  greater  than  that  at  FS313.  This  is  also  true  of  the  coefficient 
of  the  second  longitudinal-structural  mode.  The  coefficient  of  the 
third  longitudinal-structural  mode  at  FS313  is  estimated  from  data 
given  in  References  2 and  3,  since  the  generalized  shape  of  this 
mode  is  not  given.  The  coefficient  of  the  third  longitudinal-structural 
mode  has  the  same  sign  at  both  accelerometer  locations;  however,  its 
magnitude  is  four  times  greater  near  the  eg. 
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Three  different  lateral-accelerometer  locations  exist  for  the 
three  sensor  locations;  FS77,  FS186,  and  FS313.  This  sensor  is 
affected  only  by  the  first  fuselage  lateral-bending  mode,  rig:  its 
coefficient  is  negative  when  located  forward,  and  positive  when  located 
near  the  eg.  The  longitudinal-accelerometer  is  not  affected  by  the 
structural  modes  (as  discussed  previously) . 

The  pitch-rate  gyro  is  located  differently  for  each  of  the  three 
sensor  configurations  investigated;  FS77,  FS313,  or  FS383.  The  first 
longitudinal-directional  mode  coefficient  is  negative  or  equal  to  zero, 
but  becomes  more  positive  for  distances  further  from  the  nose  of  the 
aircraft.  A negative  slope  with  increasing  aft  distances  characterizes 
the  second  mode  coefficient,  and  with  a sign  change  occurring  between 
FS77  and  FS313.  The  slope  of  the  third  mode  coefficient  is  positive 
with  increasing  aft  distances,  and  with  a sign  change  at  FS313.  This 
coefficient  is  estimated  for  FS77,  due  to  a lack  of  data  in  Refer- 
ences 2 and  3. 

All  three  lateral-directional  modes  affect  the  roll-rate  gyro, 
which  is  either  at  FS77  or  FS313.  No  sign  changes  are  encountered  among 
the  mode  coefficients,  but  in  each  case  the  magnitude  decreases  with 
increasing  aft  distances.  Only  the  first  fuselage  lateral-bending  mode, 
rig,  affects  the  yaw-rate  gyro,  and  the  coefficient  of  its  effect  becomes 
more  positive  with  increasing  aft  distances.  A change  in  the  sign  of 
this  variable  occurs  between  FS313  and  FS383. 

The  parameters  selected  for  the  roll-rate  gyro  and  lateral- 
accelerometer  navigation  models  are  given  in  Tables  4-5  and  4-6.  Refer 
to  Figures  3-8  and  3-9  to  ascertain  where  these  parameters  apply  in  the 
sensor  models.  The  time  constant  of  the  lag  in  the  accelerometer  model, 
T,  -s  selected  to  be  very  small,  (10  , to  essentially  eliminate 

the  lag  from  the  sensor  model.  The  parameters  for  all  of  the  naviga- 
tion gyros  and  accelerometers  are  the  same  as  those  listed  in  Tables  4-5 
and  4-6  for  the  computer  runs  which  demonstrate  the  MIRA  simulation 
capability. 
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Table  4-5.  Ring-Laser  roll-rate-gyro  nominal  parameters. 


Parameter 

Symbol 

Value 

Misalignment  Coefficients 

Pi,  P2 

5(10“^)  rad 

Scale  Factor 

SF 

Pg 

131,328.387  pulses/rad 

Bias 

Ab 

3.006367  At  pulses  [O.OlVhl 

Pg 

Scale  Factor  Error 

e 

Pg 

' 

5(10  rad/pulse 

Table  4-6.  Quartz-flexure  lateral-accelerometer  nominal  parameters. 


Parameter 

Symbol 

1 

Value 

Misalignment  Coefficients 

’1'  ^2 

5(10"^)  rad 

Time  Constant  of  Lag 

t 

1(10“^)  s 

Scale  Factor 

SF 

n 

y 

1000  pulses/ft/s 

1 

Bias 

Ab 

0.0016  At  ft/s  [50  pn] 

Crosscoupling  Coefficients 

Q 

~ 1/ft/s  19. 6(10“^)" 

yx 

Nonlinearity  Coefficient 

C 

YY 

1.25(10"^)  1/ft/s 

Scale-Factor  Error 

e 

n 

y 

0.0001  ft/s/pulse 

4.3  MIRA  Simulation  Results 


Figures  4-2,  4-3,  and  4-4  contain  the  first  set  of  results  which 
demonstrates  the  MIRA  simulation  capability.  In  each  case,  the  sensors 
are  located  at  the  stations  recommended  for  the  SFCS  program,  as  indi- 
cated in  Table  4-1.  The  figures  differ  in  that  the  results  shown  in 
Figure  4-2  assume  perfect  sensors  with  the  vehicle  flying  in  a nontur- 
bulent  environment,  while  those  presented  in  Figure  4-3  are  obtained 
using  realistic  models  for  the  navigation  and  flight-control  sensors 
with  the  vehicle  flying  in  a lightly  turbulent  environment.  Figure  4-4 
shows  the  results  obtained  when  common  sensors  are  used  for  both  navi- 
gation and  flight  control.  Again  the  vehicle  is  flying  in  light 
turbulence. 

* This  value  is  about  10  times  larger  than  that  of  state-of-the-art 
equipment . 
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0.001  rad  ■ ; TO  ^ 0 45(10' 


Figure  4-3.  MIRA  simulation  responses — sensors  at  SFCS 
locations  (turbulent  environment) (cont. ) . 


0.0045  h/s 


Figure  4-4,  MIRA  simulation  responses — common  navigation  and  flight- 
control  sensors  at  SFCS  locations  (turbulent  environment 


The  basic  vehicle  responses  for  the  inputs  of  Figure  4-1,  are 
shown  in  Figure  4-2.  The  vehicle  achieves  a bank  angle  of  approximately 
32.5  degrees  due  to  the  application  of  the  positive  roll-stick  pulse. 

This  variable  bleeds  off  slightly  from  this  value  until  the  negative 
lateral-stick  pulse  is  applied  after  which  it  becomes  slightly  negative. 
The  yaw  angle  gradually  builds  up  to  about  37.5  degrees,  and  remains 
essentially  constant  after  the  negative  lateral  stick  is  applied.  The 
pitch  angle  decreases  slightly  from  its  trim  value  due  to  a loss  of 
lift.  Tlie  longitudinal-stick  input  is  not  large  enough  to  maintain  a 
constant  0 during  the  turning  portion  of  flight.  However,  0 increases 
substantially  with  the  application  of  the  longitudinal-stick  step  input. 

A gradual  buildup  in  velocity  occurs  as  the  pitch  angle  decreases  until 
the  longitudinal-stick  step  input  is  applied,  at  which  time  this  vari- 
able decreases  slightly.  The  angle  of  attack  is  essentially  constant 
at  its  trim  value,  until  the  longitudinal-stick  step  input,  when  it 
undergoes  a step  change.  During  most  of  the  simulation  run,  B is 
approximately  zero,  and  achieves  its  largest  values  when  the  lateral- 
pulse  inputs  are  applied.  A decrease  in  altitude  and  negative  rate  of 
altitude  change  results  during  the  turning  portion  of  flight.  The 
application  of  the  longitudinal-stick  step  input  causes  a positive  h, 
and  a subsequent  slight  increase  in  altitude  at  the  end  of  the  run. 

The  latitude,  longitude,  and  wander  angle  undergo  positive  ramp  increases 
with  time.  The  north-pointing  earth  velocity  of  the  vehicle  decreases 
slightly  from  its  trim  value,  while  the  east  component  of  the  velocity 
of  the  vehicle,  with  respect  to  the  earth,  increases  with  time  during 
the  turning  portion  of  flight,  and  then  levels  off. 

The  control-surface  deflections  and  flight-control-system  sensor 
signals  are  also  presented  in  Figure  4-2.  A pulse-aileron-surface 
deflection  .-'.id  a corresponding  pulse-roll-rate-gyro  sensor  signal  are 
produced  by  the  lateral-stick  input.  A well-damped  roll  mode  is  evident 
from  these  responses.  The  roll-to-yaw  crossfeed  and  feedback  system 
provides  a rudder-surface  deflection  which,  in  turn,  results  in  an 
essentially  constant  yaw  rate  during  the  initial  portion  of  the  time 
history.  The  initial  negative  rudder  deflection  causes  an  initial 
negative  lateral-acceleration  transient.  After  this,  the  rudder-surface 
deflection  reflects  the  envelope  of  the  lateral-accelerometer  signal, 
until  the  negative  lateral-stick  input,  when  the  time  history  repeats. 
From  these  responses,  the  Dutch  Roll  mode  appears  to  be  well-damped. 

The  effects  of  the  structural-mode  signals  are  evident  in  the  n 

y ppc 

signal  to  a small  extent.  Longitudinally,  a constant  pitch  rate 
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and  a constant  upward  normal  acceleration  are  produced  by  the 
lonyitudinal-stick  command  during  the  initial  portion  of  the  time 
history.  The  application  of  the  large  longitudinal-stick  step  input 
then  excites  the  short-period  mode.  The  high  damping  of  the  short- 
period  mode  is  evident  from  these  responses. 

The  input  variables  to  the  navigation-system  model  are  considered 
next.  Since  perfect-sensor  models  are  used  in  this  case,  these  signals 
are  repeats  of  the  flight-control-system  sensor  time  histories  with  the 
appropriate  axes  transformation  between  the  vehicle  equations  of  motion 
and  the  navigation-system  model.  The  errors  in  the  Euler  angles  are 
presented  in  Figure  4-2.  The  bank-angle  error  increases  linearly  during 
the  turning  portion  of  flight,  and  then  levels  off  near  the  end  of  thi 
simulation  run.  On  the  other  hand,  the  yaw-angle  error  remains  essen- 
tially constant  during  the  time  when  the  vehicle  is  turning,  and  then 
becomes  approximately  zero  when  the  wings  are  leveled  off.  The  pitch- 
angle  error  goes  negative  initially,  and  increases  linearly  with  time 
until  the  longitudinal-stick  input  is  removed.  After  this,  it  jumps  to 
a constant  value.  The  errors  in  longitude,  latitude,  and  altitude  all 
increase  with  time.  The  initial  offset  in  the  latitude  error  is  due  to 
the  fact  that  the  initial  value  of  the  navigation-system  latitude  is 
obtained  by  inverse  trigonometric  functions;  thus,  truncation  and  round- 
off effects  are  present.  The  error  in  v/ander  angle  is  not  presented  in 
Figure  4-2,  because  it  is  of  such  a small  magnitude.  A parabolic  shape 
characterizes  the  error  in  the  north  component  of  the  earth-relative 
velocity,  while  the  east  component  of  the  earth-relative  velocity 
increases  with  time.  The  effects  of  the  low  frequency  of  the  navigation- 
system  computations  is  evident  in  the  step-like  responses  present  in 
the  error  in  h time  history. 

A set  of  time  histories,  corresponding  to  those  of  Figure  4-2,  is 
presented  in  Figure  4-3  for  the  case  where  the  vehicle  is  flying  in 
lightly  turbulent  air  and  the  navigation  and  flight-control  sensors  are 
not  perfect.  Furthermore,  the  variables  computed  by  the  navigation- 
system  model  are  presented  in  lieu  of  the  corresponding  vehicle 
variables . 

The  major  differences  between  the  responses  of  Figures  4-2  and  4-3 
are  evident  in  the  navigation-system  input  variables  which  show  the 
effect  of  the  quantization  occurring  in  the  sensor  models  and  the  tur- 
bulence. The  turbulence  input  signals  are  also  presented  in  Figure  4-3. 

A measure  of  their  impact  on  the  sensor  signals  can  be  obtained  by 
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comparing  the  flight-control-system  sensor  responses  of  Figures  4-2 
and  4-3.  Other  time  histories,  greatly  affected  by  the  turbulence, 
are  the  sideslip  angle  and  the  rudder  deflection.  The  navigation- 
system-model  output  variables  are  presented  in  Figure  4-3  in  lieu  of 
the  corresponding  vehicle  signals.  These  variables,  as  well  as  the 
MIRA  simulation-error  variables,  closely  resemble  the  corresponding 
variables  of  Figure  4-2.  Perhaps  the  major  differences  are  evident 
in  the  h signal,  where  the  effect  of  the  low  navigation-system  compu- 
tation frequency  manifests  itself,  and  in  the  pitch-attitude  error  sig- 
nal which  is  corrupted  by  turbulence.  The  error  signals  in  Figure  4-3 
are  either  of  the  same  or  smaller  magnitude  than  the  corresponding  sig- 
nals of  Figure  4-2. 

Additional  variables  reflecting  the  errors  between  the  navigation 
perfect-sensor  signals  and  their  realistic  model  counterparts  are  also 
presented  in  Figure  4-3  to  give  the  reader  an  idea  as  to  the  magnitude 
and  shape  of  these  differences. 

Figure  4-4  contains  the  responses  obtained  when  the  navigation 
sensors  are  used  for  both  navigation  and  flight-control  functions.  The 
sensors  are  located  at  the  positions  determined  for  the  SFCS  program; 
light  turbulence  is  also  present  here.  The  most  striking  aspect  of 
Figure  4-4  is  evident  in  the  flight-control-system  signals,  which  are 
now  the  same  as  the  navigation-sensor  signals,  except  for  the  appropri- 
ate coordinate  transformation.  The  f liglit-control  sensor  signals  now 
reflect  the  effects  of  the  sensor-model  errors  and,  in  particular,  the 
quantization  which  takes  place.  However,  the  control-surface  deflec- 
tions do  not  noticeably  reflect  the  quantization  effects  present  in  the 
flight-control  signals  which  drive  the  surfaces.  The  flight-control 
system  filters  these  high-frequency  effects  out.  Once  again,  the  navi- 
gation error  signals  have  the  same  characteristic  shapes  and  magnitudes 
previously  evident  in  Figure  4-3.  The  use  of  common  sensors  (Figure  4-4) 
resulted  in  slightly  smaller  errors  in  latitude,  altitude,  and  X when 
compared  to  the  responses  of  Figure  4-3  which  were  obtained  with  inde- 
pendent sensors.  Conversely,  smaller  errors  in  <p , 6,  longitude,  and  Y 
are  achieved  using  independent  sensors.  The  errors  in  ij)  and  h are 
about  the  same  in  both  cases.  These  results  reinforce  the  feasibility 
of  the  MIRA  concept  when  considered  from  a navigation-error  point  of 
view,  although  only  limited  conclusions  can  be  drawn  due  to  the  short 
duration  of  tr.i-  simulation  runs. 
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structural-mode  excitation  is  considered  in  Figures  4-5,  4-6, 
and  4-7.  For  these  cases,  the  sensors  are  all  located  at  FS77, 

Figure  4-5  resulted  from  a nonturbulent  environment,  and  the  use 
of  perfect-sensor  models.  A turbulent  environment  and  nonperfect- 
sensor models  led  to  the  results  presented  in  Figure  4-6.  The  results 
shown  in  Figure  4-7  are  obtained  with  navigation  sensors  used  for  both 
navigation  and  flight-control  purposes  with  the  vehicle  flying  in  gusty 
air . 

In  Figure  4-5,  the  structural  modes  are  excited  only  by  the 
control-surface  deflections  with  the  lowest  frequency  longitudinal-  and 
lateral-directional  modes  excited  most.  The  generalized  bending-mode 
coefficients  (i.e.,  ...)  take  on  the  shape  of  the  dominant 

control-surface  deflection  exciting  them.  For  example,  all  of  the 
longitudinal-bending-mode  coefficients  have  the  shape  of  the  incremental 
horizontal-stabilator  deflection.  The  fuselage  first  torsion  mode,  n^, 
has  the  shape  of  the  aileron  surface  deflection,  while  the  wing  first 
asymmetric-bending  mode,  n^,  and  the  fuselage  first  lateral-bending 
mode,  n^,  and  the  fuselage  first  lateral-bending  mode,  Ug,  both  reflect 
the  shape  of  the  rudder-surface  deflection. 

The  time  histories  of  Figure  4-6  correspond  to  those  already  pre- 
sented in  Figure  4-5.  A more  constant  excitation  of  the  structural 
modes  occurs  because  of  the  turbulent  environment  in  which  the  vehicle 
is  flying.  In  the  majority  of  cases,  the  excitation  due  to  the  tur- 
bulence is  much  less  than  that  due  to  the  control-surface  deflections. 
However,  this  is  not  true  in  the  case  of  n”g  and  n^,  where  the  excita- 
tion due  to  turbulence  nearly  exceeds  that  due  to  the  control-surface 
deflections . 

Figure  4-7  contains  the  time  histories  for  the  third  member  of 
this  group  of  results  which  were  obtained  using  common  sensors  for  both 
the  navigation  and  flight-control  systems.  In  general,  the  results  of 
Figure  4-7  are  quite  similar  to  those  of  Figure  4-6.  The  most  notice- 
able difference  is  an  increase  in  the  level  of  the  structural-mode 
excitation,  which  is  particularly  evident  in  the  mode  rates  and  accel- 
erations. This  is  believed  to  be  dile  to  quantization  occurring  in 
the  navigation-sensor  models.  The  structural  modes  are  excited  by  the 
derivatives  of  the  control-surface  deflections,  and  the  effects  of  the 
quantization  are  more  pronounced  in  these  signals  than  in  the  control- 
surface  deflections  themselves.  This  result  could  have  a significant 
effect  on  the  feasibility  of  the  MIRA  concept  for  vehicles  which  tend 
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Figure  4-5.  Structural-mode  responses — ^perfect  sensors 
at  FS77  (no  turbulence) . 
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Figure  4-6.  Structural-mode  responses — sensors  at  FS77 
(turblent  environent) . 
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to  be  more  flexible  since  the  structural-mode  rates  and  accelerations 
directly  affect  the  sensor  signals. 

The  last  set  of  data  included  herein,  and  which  demonstrates  the 
capability  of  the  MIRft  simulation,  is  presented  in  Figures  4-8,  4-9, 
and  4-10.  All  of  these  responses  have  been  obtained  with  the  vehicle 
flying  in  a turbulent  environment.  The  responses  of  Figure  4-8  provide 
a basis  for  the  present  discussion,  and  have  been  obtained  with  the 
flight-control  sensors  located  at  the  SFCS  locations  and  the  navigation 
sensors  located  at  FS313,  just  forward  of  the  vehicle  eg.  For 
Figure  4-9,  common  navigation  and  flight-control  sensors  located  at 
FS313  have  been  used,  while  for  Figure  4-10,  the  common  navigation  and 
flight-control  sensors  are  located  at  FS77.  The  responses  of  Figure  4-4, 
obtained  with  common  navigation  and  flight-control  sensors  located  at 
the  SFCS  stations,  are  also  pertinent  to  the  present  discussion. 

One  difference  among  the  responses  of  the  four  cases  just  men- 
tioned is  in  the  rudder-control-surface  deflection.  The  initial  peak 
of  this  signal  becomes  less  sharply  defined  as  the  sensor  is  located 
farther  forward  of  the  eg.  The  lateral-acceleromete-'  ..  il  is  also 
different  among  these  cases  in  that  the  initial  por'  the  tran- 

sient has  different  signs.  A difference  also  exists  rn  _.ie  yaw-rate- 
gyro  output.  These  differences  arise  because  of  r and  p contributions 
to  the  lateral  acceleration,  which  are  measured  by  sensors  at  locations 
different  from  the  vehicle  eg.  The  lateral-acceleration  signal  drives 
the  rudder  surface  through  the  augmentation  system  which,  in  turn, 
affects  the  yaw-rate  response.  Certainly,  the  lateral-directional 
handling  qualities  will  be  affected  by  a change  in  the  location  of 
the  lateral  accelerometer.  However,  this  subject  has  not  been  ex- 
plored during  the  course  of  this  effort.  It  should  be  pointed  out 
that  the  designer  of  the  flight-control  system  for  a MIRA  configured 
vehicle  also  has  the  liberty  of  adjusting  the  roll-to-yaw  crossfeed 
network  of  the  augmentation  system  in  order  to  help  improve  the  hand- 
ling qualities  of  the  vehicle. 

A comparison  of  the  6^,  q,  and  An^  responses  of  Figures  4-8 
and  4-9  (for  which  the  normal  accelerometer  is  located  at  FS313)  with 
those  of  Figure  4-4  and  4-10  (for  which  the  normal  accelerometer  is 
located  at  FS77)  reveals  that  the  responses  of  the  system  with  the  more 
aft-located  normal  accelerometer  are  more  lightly  damped  than  the  others. 
The  difference  in  damping  is  due  to  the  q contribution  to  the  normal 
acceleration,  which  is  sensed  by  an  accelerometer  located  in  a position 
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Figure  4-8.  MIRA  simulation  responses — ^navigation  sensors  at  FS313 
flight-control  sensors  at  SFCS  locations  (turbulent  en 


Figure  4-9.  MIRA  simulation  responses— common  navigation  and  flight- 
control  sensors  FS313  (turbulent  environment) . 


forward  of  the  aircraft  eg.  From  the  point  of  view  of  handling  quali- 
ties, it  is  desirable  to  locate  the  normal  accelerometer  forward  of  the 
eg  in  order  to  take  advantage  of  the  increase  in  short-period  damping 
which  results.  The  fact  that  the  normal-accelerometer  location  selected 
for  the  SFCS  configuration  is  far  forward  of  the  eg  reinforces  the  idea 
that  this  is  the  best  general  location  for  this  sensor. 

A noticeable  difference  is  not  evident  among  the  roll-rate 
response  of  the  vehicle  for  the  four  cases  under  consideration.  This 
is  undoubtedly  due  to  the  fact  that  the  roll-rate-gyro  signal  is 
not  affected  by  location  to  the  extent  that  the  accelerometer  signals 
are.  The  effects  of  location  manifest  themselves  in  the  roll-rate-gyro 
signal  only  through  the  structural-mode  contributions,  which  are  nearly 
negligible  for  the  model  vehicle. 

Although  it  is  dangerous  to  draw  hard  and  fast  conclusions  from 
the  limited  amount  of  data  presented  herein,  a comparison  of  the  errors 
obtained  with  the  four  sensor  configurations  under  consideration  is 
made  in  an  attempt  to  reach  some  meaningful  conclusions.  A simple 
criterion  is  used  in  this  investigation.  The  magnitudes  of  the  maximum 
errors  and  errors  at  the  final  time  are  compared.  Equal  weight  is  given 
to  each  parameter.  Tables  4-7  and  4-8  contain  the  pertinent  error  data. 

Table  4-7.  Magnitude  of  maximum  errors  of  MIRA  simulation  runs. 


Sensor  Configuration 

Fit.  Ctl.  Sens.  Loc. 

SFCS 

FS313 

FS77 

SFCS 

Nav.  Sens. 

Loc . 

FS313 

FS313 

FS77 

SFCS 

Error 

Units 

NCS* 

CS* 

CS* 

CS* 

fi 

ft/s 

3.1215 

3.0161 

2.7541 

3.1446 

y 

ft/s 

11.4163 

11.7031 

10.2139 

10.6445 

• 

X 

ft/s 

1.1621 

1.9081 

1.1576 

1.0860 

0 

deg 

1.0013 

0.9585 

0.8815 

0.9411 

e 

deg 

0.3370 

0.4799 

0.2666 

0.3078 

deg 

0.7475 

0.7805 

0.8000 

0.7744 

\ 

deg 

5.905(10-5) 

1.054  (10-4) 

6.427  (10-3) 

6.155(10-5) 

deg 

3.942 (10-^) 

3.956(10-4) 

3.990  (10-4) 

3.896(10-4) 

h 

ft. 

31.6284 

30.6637 

23.7087 

27.2142 

* 

CS  = Common  sensors  for  navigation  and  flight  control. 

NCS  = No  common  sensors  for  navigation  and  flight  control. 
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Figure  4-8.  Magnitude  of  final  errors  MIRA  simulation  runs. 


Sensor  Configuration 

1 Fit.  Ctl.  Sens.  Loc. 

L 

SFCS 

FS313 

FS77 

SFCS 

1 

Nav.  Sens . 

Loc . 

FS313 

FS313 

FS77 

SFCS 

Error 

Units 

NCS* 

CS* 

CS* 

CS* 

I 

ft/s 

2.1651 

2.0373 

0,8269 

1.1269 

i ^ 

ft/s 

11.4163 

11.7031 

10.2139 

10.6445 

ft/s 

0.7339 

1.1114 

0,1489 

0.2453 

i 

deg 

0.9765 

0.9309 

0.8388 

0.8957 

1 ® 

deg 

0.2145 

0.4336 

0.2191 

0.2470 

1 

deg 

0.2171 

0.3436 

0.2843 

0.2249 

1 ^ 

deg 

5.660(10-5) 

1.054  (10-4) 

6.386 (10-5) 

6,087 (10-5) 

1 ^ 

deg 

3.942 (10-4) 

3.956(10-4) 

3.991 (10-4) 

3.896  (10-4) 

h 

ft 

31.6284 

30.6637 

23.7087 

27.2142 

* 

CS  = Common  sensors  for  navigation  and  flight  control. 

NCS  = No  common  sensors  for  navigation  and  flight  control. 


An  investigation  of  the  data  contained  in  these  tables  indicates 
that  the  navigation  errors  are  smallest  for  the  cases  where  the  sensors 
are  located  at  FS77  and  the  SFCS  stations.  For  example,  the  errors  in 
h and  h are  significantly  smaller  for  the  case  with  the  sensors  colo- 
cated at  FS77.  This  configuration  also  resulted  in  the  smallest  errors 
in  4'  and  X for  the  four  cases  being  considered.  These  advantages  are 
counterbalanced  because  longitude,  latitude,  and  yaw-angle  errors  are 
the  worst,  or  close  to  it„  On  the  other  hand,  the  errors  obtained  for 
the  SFCS  configured  vehicle  are  the  most  consistent  of  those  obtained, 
being  neither  the  best  nor  the  worst,  in  most  cases.  These  advantages 
are  obtained  for  the  cases  with  the  sensors  located  at  FS77  and  the  SFCS 
stations  are  so  close  in  nature  that  it  is  quite  possibly  due  to  the 
location  of  the  normal  accelerometer  at  FS77  in  both  cases,  and  the 
location  of  the  lateral  accelerometer  close  to  FS77  in  the  SFCS  con- 
figuration. 

The  configuration  with  independent  flight-control  sensors  at  the 
SFCS  locations  and  the  navigation  sensors  colocated  at  the  eg  appears 
to  be  the  next  best  system  of  those  four  considered.  The  smallest 


errors  in  latitude  and  heading  angle  are  obtained  with  this  configura- 
tion. It  is  offset  by  the  fact  that  this  configuration  resulted  in  the 
worst  altitude  and  roll-angle  errors,  and  the  next  to  the  worst  errors 
in  h, ■ Y,  and  X. 

The  configuration  with  the  common-navigation  and  flight-control 
sensors  located  at  FS313  faired  worst  of  all.  The  largest  or  next  to 
the  largest  errors  are  obtained  with  this  configuration  in  most  of  the 
variables  considered  (e.g.,  the  Euler  angles,  longitude,  h,  and  Y) , 
Significantly  larger  errors  in  latitude  and  X also  resulted.  This 
sensor  arrangement  did  not  result  in  the  smallest  errors  for  any  of  the 
variables  considered. 

The  results  of  the  investigation  indicate  that  a forward  location 
is  best  for  the  accelerometers  from  both  a navigation  and  a handling 
qualities  standpoint.  It  is  better  if  the  accelerometers  are  colocated; 
however,  some  *"reedom  is  allowed  in  the  placement  of  these  sensors  to 
minimize  t tructural-mode  pickup  with  minor  or  no  degrading  effects. 

Although  en  in  this  study,  there  is  also  a tendency  to  believe 

that  t . state  long-term  navigation  accuracy  will  be  affected 

relati  — rrttle  by  the  placement  of  the  sensors  in  the  forward  por- 
tions of  the  fuselage.  That  is,  the  q,  r,  and  p terms  measured  by 
forward-placed  accelerometers  over  a long  mission  would  be  expected  to 
be  negligible.  More  flexibility  is  allowed  in  the  placement  of  the 
rate  gyros  since  the  quantity  that  they  measure  is  dependent  upon  loca- 
tion only  to  the  extent  that  the  structural  modes  affect  them. 

4.4  Conclusions  from  MIRA  Simulation  Results 

In  summary,  the  following  major  conclusions  can  be  drawn  from  the 
data  thus  far  presented. 

(1)  The  capability  of  the  MIRA  simulation  developed  by  CSDL  has 
been  demonstrated. 

(2)  The  feasibility  of  the  MIRA  concept  is  reinforced,  since 
nearly  identical  results  are  obtained  using  independent 
and  common  sensors  for  navigation  and  flight  control. 

(3)  Forward  fuselage  locations  appear  best  for  the  accelerom- 
eters from  both  a navigation  and  handling-qualities  stand- 
point based  upon  the  simulation  runs  performed  for  this  study. 
Some  freedom  is  allowed  in  locating  these  sensors  with  re- 
gard to  their  structural-mode  contamination. 
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(4)  Flexibility  is  allowed  in  the  placement  of  the  rate  gyros, 
since  the  signals  that  they  sense  are  location-dependent 
only  to  the  extent  that  the  structural  modes  affect  them. 

(5)  The  quantization  present  in  the  navigation  sensor  models 
affects  the  level  of  structural-mode  excitation.  However, 
for  a fairly  rigid  vehicle  this  effect  does  not  manifest 
itself  to  a great  extent  in  the  sensor  signals,  control- 
surface  activity,  or  basic  vehicle  responses. 


SECTION  5 


LIFE-CYCLE  COST  MODELING  AND  MFCHANIZATION  RESULTS 


5.1  Introduction 


Life-cycle  cost  (LCC)  is  defined  as  the  total  cost  of  acquiring 
a product,  establishing  the  necessary  logistics  base  from  which  to 


deploy  and  use  the  product,  and  maintaining  the 
condition  over  some  prescribed  period  of  time, 
cycle  costing,  a program  is  often  divided  into 

product  in  operable 
For  purposes  of  life 
the  following  phases: 

(1) 

Conceptual  studies. 

(2) 

Research,  development,  test,  and  evaluation  (RDT&E) . 

(3) 

Production. 

(4) 

Operation  and  maintenance. 

(5) 

Disposal . 

The  goal  of  life-cycle  costing  (i.e.,  achieving  minimum  life-cycle 
cost)  is  depicted  in  Figure  5-1,  where  costs  are  shown  as  a function  of 
mean  time  between  failure  (MTBF) . For  purposes  of  illustration,  the 
acquisition  costs  (all  costs  incurred  through  production) , and  operation 
and  maintenance  (O&M)  costs  are  shown.  The  minimum  life-cycle  cost 
occurs  at  the  crossover  between  the  acquisition-cost  curve  and  the  O&M 
cost  curve.  Thus,  there  is  an  optimum  reliability  for  a given  piece  of 
equipment — that  reliability  which  permits  the  minimum  life-cycle  cost. 

In  the  last  decade,  increasing  emphasis  has  been  placed  on  LCC 
from  several  points  of  view. 

(1)  Minimization  of  LCC  should  be  included  as  a program  goal. 

(2)  LCC  considerations  should  b^  a factor  in  design  and  manage- 
ment decisions. 

(3)  LCC  should  be  an  independent  research  and  data  project. 

(4)  LCC  should  be  included  in  the  source  selection  criteria. 
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Figure  5-1.  The  goal  of  life-cycle  cost  analysis. 

Life-cycle  costing  is  generally  accomplished  through  the  applica- 
tion of  a diverse  spectrum  of  mathematical  models  to  generate  system 

LCC  estimates  or  to  aid  in  making  the  design  and  management  decisions 

( 9 ) 

affecting  LCC.  In  a recent  USAF  study  the  following  eight  categories 
of  LCC  models  were  identified. 

(1)  Accounting  model  — aggregates  all  costs. 

(2)  Economic  analysis  model  — treats  the  time  value  of  money 
versus  the  program  schedules. 

(3)  Cost-estimating  relationship  model  — relates  cost  tc 
parameters . 

(4)  Reliability/improvement  cost  model  — relates  cost  to 
reliability. 

(5)  Level-of-repair  analysis  model  — determines  minimum-cost- 
maintenance  policy. 


92 


(6)  Maintenanc  -manpower  planning  model  — evaluates  alternative 
maintenance-manpower  requirements  or  effects  of  alternative 
equipment  designs  on  maintenance-manpower  requirements. 

(7)  Inventory  management  model  — determines  spare-part  stock 
levels  to  minimize  spares  costs. 

(8)  Warranty  model  — evaluates  government  maintenance  versus 
contractor  maintenance  under  warranty. 

5. 2 Objective 

The  goal  of  the  life-cycle-costing  effort  in  MIRA  is  the  develop- 
ment of  LCC  techniques  as  design  tools.  The  traditional  approach  in  the 
application  of  LCC  is  to  sum  all  the  costs  which  would  be  included  if  a 
proposed  equipment  design  were  put  into  production.  In  this  context, 

LCC  becomes  an  important  source  selection  criterion  and  in  some  instances 
may  be  the  most  important  criterion.  However,  most  of  the  decisions 
which  affect  LCC  have  been  made  during  the  early  phases  of  the  program; 
in  particular,  during  conceptual  studies  and  engineering  design.  There- 
fore, it  is  most  appropriate  that  part  of  the  MIRA  effort  should  address 
the  development  of  LCC  techniques  to  aid  the  system  designers  in  working 
towards  a design  which  will  meet  performance  requirements  at  an  afford- 
able cost  to  the  government,  rather  than  waiting  to  determine  the  costs 
after  a design  has  been  chosen. 

5 . 3 LCC  Model  Survey 

With  respect  to  MIRA,  the  most  useful  of  the  models  identified 
in  Section  5.1  is  the  accounting  model  which  estimates  the  aggregate 
of  all  the  costs,  and  yet  permits  preselected  cost  areas  to  be  consid- 
ered separately.  As  the  program  work  progresses,  such  an  approach 
would  be  refined  to  produce  increasingly  better  estimates  of  LCC.  How- 
ever, the  application  of  one  or  more  of  the  other  types  of  models  is 
not  ruled  out,  and  may  be  required  as  the  MIRA  effort  progresses. 

The  first  item  considered  in  the  LCC  effort  for  MIRA  was  the 

review  of  available  LCC  accounting  models  from  both  DoD  and  industry. 

(9) 

A USAF  review  of  LCC  models  developed  for  and  by  the  military  con- 

cludes that  most  of  the  general  models  (including  the  accounting  type) 
are  inadequate  because  they  lacked  focus  on  specific  decisions  or  trade- 
offs, did  not  describe  a particular  equipment  or  system  very  well,  and 
required  excessive  and  specialized  data.  An  earlier  review  of  LCC 
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models  by  industry indicated  that  the  available  models  relating  to 
military  systems  are  too  limited  in  scope,  too  complex,  and  required 
unavailable  data,  or  are  excessively  time  consuming. 

Two  comprehensive,  accounting-type  LCC  models  have  been  selected 
for  possible  application  to  MIRA.  Both  have  been  recently  developed  by 
knowledgeable  LCC  personnel  and  avoid  many  of  the  pitfalls  described  in 
the  two  surveys  referenced  previously.  They  are  as  follows; 

(1)  The  "Cost  Reduction  Is  Everyone's  Responsibility"  (CRIER) 
model,  developed  by  the  Life-Cycle-Cost  Task  Group  of  the 
Joint  Services  Data  Exchange  for  Inertial  Systems.  This 
Task  Group  is  composed  of  representatives  from  the  Armed 
Services  (Air  Force,  Navy,  and  Army),  industry,  and  the 
academic/research  community.  CRIER  is  very  comprehensive 
and  includes  R&D,  acquisition  and  operations,  and  main- 
tenance factors.  It  is  available  without  charge  and  recom- 
mended by  the  Task  Group  as  the  standard  LCC  accounting-type 
model  for  the  Armed  Services. 

(2)  The  "AGMC  Model  for  Inertial  Navigation  Systems,"  developed 
by  personnel  at  the  USAF  Aerospace  Guidance  and  Metrology 
Center.  In  an  overall  sense,  this  model  is  less  comprehen- 
sive than  CRIER,  but  also  includes  R&D,  acquisition,  and 
O&M  costs,  as  well  as  being  available  without  charge. 

5.4  The  CRIER  Model 

The  CRIER  model  was  the  only  one  considered  during  this  phase  of 
the  MIRA  effort  since  it  is  an  accounting-type  model  developed  specifi- 
cally for  inertial  systems.  It  is  also  applicable  to  all  types  of 
electronics  systems  and  comprehensive,  including  all  costs  from  concep- 
tual studies  through  operations  and  maintenance  to  disposal.  The  CRIER 
model  can  be  used  in  design  tradeoff  studies  and  is  readily  available. 
CSDL  personnel  are  familiar  with  CRIER  since  they  participated  in  the 
generation  of  the  rnodeT  through  membership  in  the  Life-Cycle-Cost  Task 
Group. 

Some  of  the  major  advantages  of  CRIER  are: 

(1)  It  includes  the  interests  of  both  industry  and  the  military. 

(2)  It  allows  for  general  or  detailed  data  inputs. 
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(3)  It  contains  general  cost  factors  which  can  be  overridden  as 
specific  data  become  available. 

(4)  It  provides  selectable  standard  output  reports. 

(5)  It  allows  the  user  to  develop  his  own  output  reports. 

(6)  It  contains  override  features  when  only  specific  sections 
of  the  model  are  of  interest. 

Special  character.-Stics  of  CRIER  permit  the  user  to  select  appli- 
cable cost  categories  and  select  the  level  of  input-data  detail.  Also, 
input  data  are  required  only  for  applicable  cost  categories,  and  spares 
prediction  is  cost  optimized  using  Poisson  distribution  probabilities. 

Table  5-1  lists  the  major  cost  elements  included  in  CRIER  and 
illustrates  the  comprehensive  nature  of  this  model. 

The  CRIER  LCC  model  has  been  installed  on  the  CSDL  AMDAHL  470 
central  computer,  and  is  in  the  process  of  checkout  and  evaluation. 

5 . 5 LCC  Impact  on  Design 

The  traditional  approach  to  life-cycle  costing  is  illustrated  in 
the  upper  half  of  Figure  5-2.  The  requirements  for  the  equipment  are 
translated  into  a set  of  detailed  performance  specifications,  and  an 
engineering-design  LCC  analysis  is  then  performed  to  obtain  the  total 
costs  of  the  proposed  design. 

In  the  MIRA  approach,  illustrated  in  the  lower  half  of  Figure  5-2, 
LCC  analysis  becomes  an  integral  part  of  the  performance-specification 
and  engineering-design  efforts  by  providing  continuous  cost  feedback 
during  these  processes. 

After  a specific  approach  to  equipment  design  is  selected,  the 
detailed  equipment  and  LCC  designs  proceed  with  continuous  feedback 
from  the  LCC  efforts  to  the  engineering  effort.  Thus,  tradeoffs  are 
conducted  during  all  phases  of  performance  specification  and  equipment 
design  (from  major  decisions,  to  lower  level  detailed  design  and  O&M 
concepts ) . 

Tradeoff  studies  should  specifically  address  elements  which 
experience  has  shown  to  be  high-cost  drivers  in  order  to  achieve  the 
optimum  design  approach  at  minimum  LCC.  Reliability  is  a major  cost 
driver,  and  affects  the  costs  of  spares,  maintenance,  transportation, 
etc.  The  CRIER  model  enables  tradeoff  studies  to  be  conducted  as  basic 
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Table  5-1.  CRIER  major  cost  elements 


Aerospace  Ground  Equipment/Ground-Support  Equipnent/Test  Equipment 


TRADITIONAL 


MIRA 


Figure  5-2.  Approaches  to  life-cycle  costs. 

cost  input  data  are  developed.  Typical  examples  of  LCC  impact  areas 
which  could  be  subjected  to  tradeoff  analysis  include; 

(1)  Performance  versus  reliability. 

(2)  Performance  versus  test/test-equipment  cost. 

(3)  Built-in  calibration  versus  calibration  plus  software  change 

(4)  Built-in  test  versus  test  by  software  versus  troubleshooting 

(5)  Throwaway  versus  repair. 

(6)  New  parts  versus  standard  parts. 

(7)  Redundancy  levels. 

(8)  Plug-in  modules  versus  hard-wired  modules. 

(9)  Reliability  versus  maintainability. 

(10)  Optimum  maintenance  concepts. 

The  integration  of  LCC  techniques  into  the  design  process  can 
significantly  aid  in  achieving  a MIRA  design  which  will  meet  performance 
requirements  at  a minimum  life-cycle  cost. 
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SUMMARY  AND  CONCLUSIONS 


6.1  Phase  I Sununary 

The  Phase  I MIRA  effort  at  CSDL  was  dedicated  toward  developing 
tools  and  techniques  which  will  be  useful  to  the  Air  Force  in  the  design 
and  evaluation  of  MIRA  systems  over  the  long  term.  It  was  felt  that  the 
multidisciplinary  goal  of  the  MIRA  program  required  the  development  of 
a unique  capability  which  could  evaluate  all  aspects  of  candidate  MIRA 
systems. 

The  primary  tool  intended  tor  the  MIRA  program  was  an  all-digital 
real-time  dynamic  and  kinematic  simulation,  which  could  be  used  to 
evaluate  a wide  class  of  MIRA  systems  in  many  different  environments. 

Although  a hybrid  (analog/digital)  simulation  would  have  repre- 
sented a more  classical  approach  to  the  problem,  particularly  since 
high-frequency  dynamics  (e.g.,  body-bending  response)  are  of  interest, 
it  was  felt  that  an  all-digital  approach  provided  a substantial  advan- 
tage in  terms  of  flexibility,  stability,  and  repeatability.  The  digital- 
processor  technology  available  at  CSDL  could  support  an  all-digital 
real-time  simulation,  and  has  significant  growth  potential  to  support 
future  program  developments. 

Sufficient  simulation  capability  was  required  to  include  at  least 
the  lower  frequency-bending  dynamics,  and  provide  a realistic  dynamic 
environment  for  the  MIRA  instruments.  Although  it  was  anticipated  that 
these  effects  would  have  a relatively  small  impact  on  the  navigation 
performance,  it  was  felt  that  they  would  affect  the  flight-control  per- 
formance significantly. 

All  of  the  simulation  goals,  except  for  the  full  real-time  capa- 
bility, were  realized  in  Phase  I.  The  lack  of  a real-time  capability 
would  represent  a major  handicap  for  a program  which  was  heavily 
oriented  toward  testing,  but  had  a relatively  small  impact  on  this 
developmental  portion  of  the  program. 
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Thirteen  different  MIRA  configurations  were  investigated;  they 
represented  a wide  range  of  system  configurations  and  constraints. 
Although  the  test  periods  used  in  this  phase  were  short,  the  systems 
demonstrated  comparable  performance  for  both  navigation  and  flight  con- 
trol in  the  tests  accomplished.  Although  strong  conclusions  were  not 
drawn  concerning  the  relative  merits  of  these  configurations,  it  was 
felt  that  these  cases  did  successfully  demonstrate  the  flexibility  of 
the  basic  tool  to  accommodate  a wide  class  of  MIRA  systems. 

The  vehicle  model  used  in  the  Phase  I simulation  was  well 
documented,  and  a significant  amount  of  equivalent  data  were  available 
to  CSDL  from  the  literature.  This  comparative  data  provided  a valuable 
check  during  the  development  phase  of  the  simulation.  Although  not  a 
true  state-of-the-art  vehicle,  the  performance  was  sufficient  to  point 
out  some  of  the  limitations  of  strapdown  systems  when  used  in  the  MIRA 
context.  Conversion  of  the  vehicle  model  to  other  aircraft  during  this 
phase  was  limited  only  by  data  availability,  particularly  for  data  re- 
lating to  body-bending  dynamics. 

The  relatively  straightforward  stability-augmentation  control 
system  appeared  to  function  adequately  with  either  dedicated  flight- 
control  sensors  or  flight-control  data  derived  from  navigation  sensors. 
In  systems  with  very  small  phase  margins,  the  derived  data  might  prove 
unacceptable,  particularly  if  phase  stabilization  of  the  body  bending 
is  required. 

The  navigation  algorithms  included  in  the  simulation  were  derived 
from  the  Standard  Navigation  Software  effort,  and  they  proved  adequate 
for  the  system  considered.  Because  of  the  relative  magnitudes  involved, 
however,  evaluation  of  navigation  errors  was  difficult,  and  sensor 
errors  were  indistinguishable  from  model,  algorithm,  and  numerical 
roundoff  effects.  A modification  to  the  navigation  algorithms  to  allow 
a clearer  evaluation  of  sensor  effects  is  proposed  for  Phase  II. 

The  life-cycle-cost  effort  provided  survey  data  concerning  life- 
cycle-cost models  and  some  preliminary  parametric  data  concerning  the 
impact  of  design  decisions  upon  life-cycle  cost.  This  overall  effort 
was  hampered  by  the  highly  empirical  nature  of  most  life-cycle-cost 
analysis  models,  and  the  high  degree  of  specificity  required  for  many 
system  parameters.  Without  a specific  system  design,  the  impact  of 
many  of  these  factors  could  not  be  assessed.  It  is  felt  that  the  over- 
all concept  of  utilizing  life-cycle-cost  data  on  an  a priori  basis  is 
promising  for  an  effort  such  as  MIRA,  and  the  communication  channels 
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between  the  engineering  and  life-cycie-cost  areas  should  be  kept  open 
to  provide  coordinated  system  design  decisions. 

In  summary,  it  is  felt  that  the  Phase  I effort  achieved  its  primary 
goals  of  demonstrating  a viable  approach  for  evaluating  MIRA-type  avion- 
ics systems  from  many  aspects.  The  basic  tools  developed  during  this 
effort  are  currently  available  at  CSDL  and  will  be  used  in  all  further 
efforts  in  this  area. 

6 . 2 Phase  II  Inputs 

During  the  course  of  the  Phase  I MIRA  effort,  there  were  several 
areas  of  investigation  encountered  which  appeared  to  suggest  meaningful 
candidates  for  the  follow-on  effort.  Some  of  these  areas  were  antici- 
pated at  the  beginning  of  the  Phase  I effort,  while  others  only  became 
apparent  as  the  effort  progressed. 

They  are  as  follows: 

(1)  Real-Time  Capability  — It  is  highly  desirable  to  provide 
full  real-time  capability  for  the  complete  simulation. 

This  not  only  allows  man-in-the-loop  activity,  but  greatly 
improves  the  efficiency  of  the  navigation  testing  which,  by 
its  very  nature,  must  represent  relatively  long  duration 
tests. 

(2)  Alternate  Aircraft  Models  — The  basic  simulation  structure 
is  adaptable  to  a wide  class  of  vehicle  models.  A more 
state-of-the-art  fighter  model  is  desirable  from  a system- 
application  viewpoint,  and  a transport  model  may  be  desir- 
able from  a MIRA  program  validation  standpoint.  The  largest 
constraint  in  adapting  the  simulation  to  different  vehicle 
models  is  obtaining  sufficient  data,  particularly  in  the 
area  of  bending  dynamics.  It  is  not  felt  that  there  would 
be  significant  gain  in  converting  to  a rigid-body  model  if 
bending  data  is  not  available. 

(3)  Alternate  Sensor  Models  — The  simulation  is  readily  adapted 
to  utilization  of  alternate  sensor  models  for  MIRA  systems. 
The  approach  of  modeling  sensor  errors  is  quite  general, 
and  there  exist  data  on  a wide  class  of  suitable  sensors. 

(4)  Navigation  Error  Analysis  — The  approach  to  navigation 
error  analysis  should  be  modified  to  that  discussed  in 
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Section  3.13  in  order  to  allow  more  efficient  separation  of 
sensor  errors  from  all  other  effects  detrimental  to  naviga- 
tion accuracy. 


(5)  Study  of  Aided  Navigators  — From  a system  design  viewpoint, 
it  appears  highly  desirable  to  expand  the  MIRA  systen:  con- 
cepts studied  to  include  a variety  of  aided  navigators 

if 

(e.g.,  TACAN  or  GPS)  . Such  aided  systems  might  have  sig- 
nificantly different  MIRA  sensor  requirements  than  unaided 
systems . 

These  items  seem  to  represent  areas  of  natural  evolutionary 
development  for  the  MIRA  effort  at  CSDL  to  allow  its  utilization  as  a 
meaningful  evaluation  tool  for  potential  MIRA  systems  and  for  testing 
and  developing  MIRA  system  concepts  of  interest  to  the  Air  Force. 


* TACAN  = Tactical  Aid  to  Navigation. 
GPS  = Global  Positioning  System. 
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